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1.  INTRODUCTION 


This  manual  Is  divided  into  two  parts.  Section  i  constituted  the  user 
section,  describing  the  overall  flow  of  the  codes,  input  data  requirements, 
and  parameter  relations.  At  the  end  of  the  section,  a  series  of  sample 
problems  is  given  to  illustrate  the  use  of  the  codes.  The  sample  problems 
exercise  all  parts  of  the  A-STAR  codes  and  serves  as  a  check  case  for  the 
installation  of  the  codes  on  a  given  machine.  The.  problems  are  arranged  ia 
increasing  computational  complexity.  Section  3  provides  c.  detailed 
description  of  all  codes  and  subroutines.  The  user  need  not  refer  to  Section 
3  or  to  the  analysis  in  Volume  I  of  this  report  to  successfully  install  and 
execute  the  programs  on  a  particular  computer  system. 


1.1  Formulation  and  Method  of  Solution 


A  detailed  description  of  the  theoretical  formulation  forming  the  basis 
of  the  A-STAR  programs  is  given  in  Volume  I.  The  analysis  is  based  on  the 
solution  of  the  generalized  electric-field  integral  equation  (EF1E)  for  all 
the  BOT  configurations  treated  with  the  method  of  moments  (MM.'  technique.  The 
unknown  currents  on  the  surface  of  the  body  and  the  antennas  are  expanded  in  a 
series  of  expansion  functions  appropriate  to  each  of  the  major  parte  of  the 
BOT  configuration:  the  BOT  surface,  the  caps,  the  wire  radiators,  the 
junction,  and  the  edge  transition  regions  on  the  body.  The  formulation  allows 
the  effect  of  some  or  all  of  these  components  to  be  treated.  The  resulting 
analysis  and  implementing  computer  programs  provide  a  unified  treatment  of  a 
broad  class  of  radiation/scattering  problems  associated  with  radiators  on 
asymmetric  surfaces. 


1.2  Program  Features 

A  partial  synopsis  of  salient  features  and  options  of  the  A-STAR  codes  is 
given  below: 

Flexible  body  geometry  -  The  algorithm  can  treat  a  finite-length  BOT  of  any 
asymmetric  cross  section.  The  BOV  can  be  open  or  closed  (i.e.,  capped).  The 
caps  are  assumed  to  be  planar  and  perpendicular  to  the  axis  of  translation  of 
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the  BOT.  Degenerate  forms  of  a  BOT  configuration  can  also  be  analysed  with 
this  algorithm  each  as  flat  plates,  parabolic  and  square  cylinders,  and  disks 
(circular  or  noncircular). 

Arbitrary  aperture  antenna  placement  and  excitation  -  The  algorithm  treats 
single  or  multiple  rectangular  aperture  (slot)  antennas  embedded  in  the  BOT 
surface.  Location  of  the  antennas  can  be.  anywhere  on  the  BOT  except  the 
caps.  Spacing  between  adjacent  apertures  can  be  electrically  small.  All 
apertures  can  be  asymmetric. 

Arbitrary  off-surface  (wire)  antenna  placement  and  excitation  -  A  wide  variety 
of  wire  antenna  configurations  such  as  monopoles  and  loops,  both  active  and 
passive,  can  be  treated  with  these  codes.  The  antennas  can  be  anywhere  on  or 
near  the  BOT  surface.  (Location  of  an  antenna  on  the  caps  is  excluded.)  The 
user  can  specify  the  use  of  a  special  junction  representation  for  the  antenna 
attachment  points  on  the  80T,  which  Js  of  special  importance  in  the  treatment 
of  parasitic  elements. 

Surface  currents  -  The  algorithm  out, uts  the  surface  currents  on  the  BOT 
surface  in  spatial  and  modal  form.  Both  the  magnitude  and  the  phase  of  the  t- 
and  z-directed  components  are  given. 

Choice  of  polarization  and  radiation  planes  -  The  radiated  fields  and  the 
power  gain,  normalized  to  an  isotropic  radiator,  are  given  in  the  user- 
specified  plan  for  6  and  polarization. 

Arbitrary  choice  of  sampling  points  for  near  fields  -  The  electric  and 
magnetic  fields  (six  components  in  all)  resulting  from  currents  induced  by 
antennas  on  a  general  BOT  or  induced  by  incident,  fields  are  computed  at  user- 
specified  points  in  the  vicinity  of  the  BOT  surface. 

User-oriented  program  features  -  The  programs  are  liberally  interspersed  with 
comment  statements  referring  to  Volumes  I  and  II.  The  codes  are  modular,  and 
error  checking  is  provided  for  the  user  at  critical  steps  in  the  computation 
sequence.  The  input  data  requirements  are  minlmired.  All  inputs  are  printed 
out  for  user  verification. 
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Sample  problems  -  The  unique  features  of  the  algorithm  are  demonstrated  with 
sample  problems  in  Section  2.6.  The  problems  have'  been  tested  and  compared 
with  corresponding  results  using  other  methods  and  available  experimental 
data. 

As  in  all  mathematical  and  numerical  modeling,  caution  must  be  exercised 
in  applying  the  analysis  to  certain  cases  where  there  are  abrupt 
discontinuities  in  the  BOT  surface  and  at  resonance  conditions.  To  achieve 
sufficient  accuracy,  a  large  number  of  modes  and  surface  (triangle)  expansion 
functions  may  be  necessary  for  these  problems. 


2.  USSR  SECTION 


The  A-STAR  computer  ..odes  are  written  In  FORTRAN  IV,  consisting  of  e^x 
impedance  generation  programs  (i.e.,  BOTZSS,  BOTZSW,  BOTZWW,  BOTZSC,  BOTZCC, 
and  BOTZCW) ,  two  matrix  Inversion  programs  (i.e.,  BOTINV  and  BOTINVA),  four  EH 
signature  prediction  codes  (i.e.,  EOTAC,  BOTRA,  BOTSCM,  and  BOTSCB),  and  one 
utility  program  (BOTSEG)  as  shown  in  Figure  1.  All  of  the  programs  have  user- 
oriented  inputs,  which  can  be  generated  easily  for  a  given  problem.  In 
addition,  binary  files  are  generated  for  the  impedance  matrices  and  the 
inverted  matrices  by  several  of  the  programs.  Examples  demonstrating  the 
input  data  and  output  from  each  of  the  programs  are  discussed  in  the  sample 
problems  (Sec'-lon  2.6). 

2.1  Implementation  of  Computer  Codes 

All  of  the  A-STAR  computer  codes  are  writcen  in  USA  Standard  FORTRAN  IV, 
with  the  exception  of  end-of-file  (EOF)  checks  in  programs  BOTH."/,  BOTINVA, 
BOTAC,  BOTRA,  BOTSCM,  and  BOTSCB.  The  Eue  checks  given  in  the  code  listings 

in  Volume  III  are  specific  to  the  compiler  used  in  the  program  development 

(i.e.,  CDC  CYBER  175  system).  The  functioning  of  the  EOF  checks  in  the 
present  listings  is  as  follows; 

IF  EOP(u)  a,  b 
u  -  unit  number  to  check. 

a  -  statement  label  to  branch  to  if  an  EOF  is  encountered, 

b  -  statement  label  to  branch  to  if  an  EOF  is  not  encountered. 

The  device  unit  numbering  convention  is  as  follows: 
unit  number  5  -  card  reader 
unit  number  6  -  line  printer 

unit  number  7  -  binary  output  file.  Each  of  the  six  impedance-generating 
.  .. -es  creates  a  new  binary  output  file  on  unit  number  7  for  storing 
impedance  matrices.  The  two  inversion  codes  create  a  new  binary  output 
file  on  unit  number  7  for  storing  the  new  inverted  system  matrix  (see 
Figure  1), 

unit  number  1  -  binary  input  file.  Each  of  the  four  EM  signature  codes 
reads  the  inverted  system  matrix  from  unit  number  1.  The  BOTINVA 


inversion  code  reads  the  old  inverted  system  matrix  from  unit  number 
1.  The  BOTINV  inversion  code  reads  the  BOT-BOT  impedance  matrix  f*-om 
unit  number  1  (see  Figure  1). 

unit  number  2  -  binary  input  file,  which  is  read  by  the  BOTINVA  code.  It 
should  contain  the  BOT-wire  impedance  matrix  when  wires  are  being  added 
to  the  old  system  matrix,  and  the  BOT-cap  impedance  matrix  when  caps 
are  being  added. 

unit  number  3  -  binary  input  file,  which  is  read  by  the  BOTINVA  code. 

It  should  contain  the  cap-wire  impedance  matrix  if  caps  and  wires  are 
contained  In  the  new  system  matrix. 

unit  number  4  -  binary  input  file,  which  is  read  by  the  BOTINVA  code.  It 
should  contain  the  wire-wire  impedance  matrix  when  wires  are  being 
added  to  the  old  system  matrix,  and  the  cap-cap  Impedance  m  trix  when 
caps  are  being  added. 


U»er  Specified  Input* 


Usar  Specif  ied  Outputs 


Antenna  coupling  analytic  BOTAC 


«  Wire  an'.ennas 
•  Slot  antenna* 


Calculation  of  radiation  patterns:  BOTRA 

•  Principal  polarization/planes 

•  Slot  antennal 

•  Wire  antennas 

•  Current  distribution  on  body 

•  Near-field  analysis 


Calculation  of  bistatic  scattering:  BOTSCB 

•  Vertical/horizontal  polarization 

•  Arbitrary  scattering  angle* 

•  Current  distribution  on  body 

•  Near-field  analysis 


Calculation  of  monostatic  scattering:  BOTSCM 


•  Vertical/horizontal  polarization 

•  Arbitrary  scattering  angles 


1.2  BOTSF.G  Descrlpticn 

BOTSEG  is  a  utility  program,  which  can  be  used  to  segment  the  BOT 
generating  curve  using  a  limited  set  of  data  points  (see  coordinate  geometry 
in  Figure  Is,.  The  required  input  data  and  formats  are  described  below. 
READ(5,1)NPT3,NP 

1  F0RMAT(2I3) 

NPTS  -  Number  of  input  data  points  used  to  describe  the  BOT 
generating  curve. 

NP  -  Number  of  equally  spaced  BOT  gene rating-curve  data  points  to  be 
calculated.  (The  calculated  data  are  used  as  input  to  the  remaining 
BOT  programs.) 

DO  100  I  «  1,  NPTS 

100  REAT( 5, 2)XTAB( l) , YTAB(I) ,XC( I) , YC( I) 

2  FORMAT(4E10.4) 

XTAB(I)  -  x  coordinate  of  the  I-th  point  on  the  input  curve  (meters). 
YTAB(I)  -  y  coordinate  of  the  1-th  point  on  the  input  curve  (meters). 
XC(I),YC(I)  -  (meters)  indicate  whether  the  points  [XTAB(I),  YTAB(I)] 
and  [XTAB(I+l),YTAB(J+i)J  are  connected  by  a  straight-line  segment  or 
an  arc  with  changing  radius.  If  XC(1)  -  YC(I)  -  0,  the  segment  is 
straight.  Otherwise,  [XC(I) ,YC(I) J  is  assumed  to  be  the  center  of  an 
arc  with  the  above  end-points,  where  the  radius  changes  linearly  with 
angle  from  [ XTAB ( I ) , YTAB ( I ) ]  to  [XTAB(I+1) ,YTAB(I+1) ] ,  subtending  the 
angle  of  the  triangle  formed  by  the  three  points.  If  the  three  points 
are  colllnear,  the  direction  of  the  arc  is  clockwise. 

Example :  Consider  a  cylinder  with  a  cross-section  given  in  Figure  2.  The 
generating  curve  for  this  body  can  be  represented  with  three  points  (NPTS  “  3) 
as  follows: 

0  0  0  0 

3  13  0 

4  0  0  0 

BOTSEG  is  currently  dimensioned  to  handle  100  points  on  the  input  curve 
and  83  points  on  the  output  curve.  The  NP  data  points  used  to  represent  the 
BOT  generating  curve  determine  the  location  and  number  of  triangle  functions 
and  their  derivatives  (denoted  by  arrays  T  and  TP  in  the  codes)  used  to 
discretize  the  unknown  currents  on  the  BOT  surface.  As  an  example,  Figure  3 
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shows  a  portion  of  a  BOT  generating  curve  defined  by  points  t^,  •  . 

t1  with  the  triangle  functions  centered  at  t^,  t^,  etc.  The  triangle 
functions  span  five  data  points,  with  adjacent  functions  overlapping. 


Input  Data  Description  and  Formats 


A  detailed  description  of  the  input  data  required  by  a.’l  A-STAB.  codes 
(except  the  inversion  codes,  as  described  in  Section  2.3.8)  is  given  below, 
READ  statements  and  formats  are  listed  in  the  order  in  which  they  appear 
within  the  program,  followed  by  a  description  of  the  required  data. 


REAP(5, 1)BK 
1  FORKAT(E15.7) 

BK  -  Wave  number  for  the  problem  (meters-*). 


Figure  2.  Repretentation  of  «  BOT  with  a  cron  taction 
formed  by  a  wedge  and  arc  augment. 


Figure  3.  Triangle  function*  on  BOT  turf  ace. 


2.3.1  BOT  Data  Set 


This  data  set  is  required  by  all  programs. 

READ(5,2)  NMODE, NSP, NBAND 

2  FORMATC3I3) 

NMODE  -  Number  of  nonnegative  modes  to  be  considered  (i.e.,  there  will 
be  2*NM0DE  -  1  total  modes). 

NSP  -  Number  of  diagonal  bands  to  be  used  in  each  ZCt  ,  Z*C  , 

m  >  n  ni  *  n 

f  2  ZZ 

Z  ,  and  Z  submatrix.  NSP“1  indicates  that  only  the 
m,n’  m,n 

diagonal  terms  are  nonzero  in  each  submatrix;  NSP»2  indicates 
that  only  diagnonal  and  off-diagonal  terms  are  nonzero,  etc. 

If  NSP  >  (NP-3)/2,  each  submatrix  is  full,  where  NP  is 
described  below. 

NBAND  -  Number  of  submatrix  diagonal  bands  to  be  calculated  by 
BOTZSS.  NBAND-i  Indicates  that  only  the  diagonal  Zm  n 
submatrices  ate  to  be  calculated.  NBAND-2  indicates  that 
diagonal  and  off-diagonal  Z  n  submatrices  are  to  be 
calculated,  etc.  If  NBAND  >  2*NM0DE-1,  the  entire  ZBQT  matrix 
is  calculated. 

READ(5,3)NP 

3  F0RMATU3) 

NP  -  Number  of  points  used  to  describe  the  BOT  generating  curve. 

(NP  must  be  odd.)  If  the  BOT  generating  curve  is  closed  (i.e., 
the  first  and  last  points  coincide),  the  programs  will  increase 
NP  by  two  and  add  two  points  to  the  generating  curve  (i.e.,  the 
YB  and  XB  arrays  described  below).  This  new  NP  should  be  used 
in  all  definitions  involving  NP  (e.g.,  dimensions). 

READ(5,4)(YB(I) ,1*1 ,NP) 

4  FORMAT( 10F8.4) 

YB  -  Array  of  y  coordinates  for  the  generating  curve  (meters). 

READ(5,5)(XB(I),I-1,NP> 

5  FORMAT( 10F8.4) 

XB  ••  Array  of  x  coordinates  for  the  generating  curve  (meters). 
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READ(5,6)BL 

6  FORMAT(F8.4) 

BL  -  Half  length  of  the  BOT  (meterB). 

2.3.2  Cap  Data  Set 

This  data  set  Is  required  by  all  programs  and  Is  used  for  putting  end- 
caps  on  the  open  BOT  (see  Figure  4). 

READ(5 ,7)NC, NPR,NE 

7  FORMAT(3I3) 

NC  -  Number  of  flat  end-caps  on  the  BOT. 

NPR  -  Number  of  radial  points  used  to  describe  each  cap.  (NPR  must 

be  odd.) 

NE  -  Indicates  whether  or  not  the  BOT/cap  transition  term  is  •"»  be 
Included.  NE-0  Indicates  no  transition  terms,  and  NE  *  0 
Indicates  that  the  BOT/cap  transition  term  should  be  used. 


Figure  4.  Representation  of  •  capped  BOT. 
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If  NE  *  0,  the  programs  set  NE"NC.  This  new  NE  should  be 
used  in  all  definitions  Involving  NE  (e.g.,  dimensions). 

If  NC=0,  che  rest  of  thic  data  set  is  ignored. 

READ(b,8)XC,YC 

8  FORMAT(2F8.4) 

XC  -  x  coordinate  of  the  cap  center  (meters). 

YC  -  y  coordinate  of  the  cap  center  (meters). 

READ(5,9)(ZC(I),I-1,NC) 

9  F0RMAT(10F8.4) 

ZC  -  Array  of  z  coordinates  for  the  caps  (meters). 

RLAD( 5 , 10) (RHOC( I ) , I-l , NPR) 

10  FORMAT( 10F8.4) 

RHOC  -  Array  of  normalized  radial  coordinates  used  to  describe 
the  caps. 

(0<  RH0C(I)<RH0C(I+iW.) 

If  NE  *0,  the  following  is  required. 

READ(5 , 1 1) (ZE( I) , T«1 ,NC) 

11  F0RMAT( 10F8.4) 

ZE  -  Array  of  z  coordinates  that  specifies  the  starting  z  location 
for  the  BOT/cap  transition  term  (meters).  (-BL<ZE(I)<BL. ) 

2.3.3  Nlre  Data  Set 

This  data  set  is  required  by  all  programs. 

READ(5, 12)NW,NPW,NJ 

12  F0RMAT(3I3) 

NW  -  Number  of  wire  antennas  on  the  BOX.  Each  antenna  must  be 
represented  by  an  odd  number  of  points. 
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NPW  -  Total  number  cf  data  points  representing  all  of  the  wire 
antennas  on  the  BOT. 

NJ  -  Number  of  BOT/wire  junction  pointB  to  be  Included. 

If  NW-0,  the  rest  of  this  data  Bet  Is  Ignored. 

READ(5,l3)(XW(I),l«l,NPW) 

13  FORMAT( 10F8.4) 

XW  -  Array  of  x  coordinates  for  the  wires  (meters). 

READ(5,14)(YW(I),I«*1,N?W) 

14  FORMAT ( 10F8.4) 

YW  -  Array  of  y  coordinates  for  the  wires  (meters). 

READ(5,15)(ZW(1),I«1,NPW) 

15  FORMAT( 10F8.4) 

ZW  -  Array  of  z  coordinates  for  the  wires  (meters). 

RZAD(5,1S)(INDW(I),1-1,NW) 

16  FORMAT( 1013) 

INDW  -  Array  containing  the  starting  index  for  each  antenna  in  the  XW, 
YW,  and  ZW  arrays.  INDW ( I)  is  the  starting  index  of  the  1-th 
antenna. 

READ(5,l7)(RADW(I),I«i,NW) 

17  F0RMAT( 10F8.4) 

RADW  -  Array  containing  the  redii  of  the  antennas.  RADW(I)  Is  the 
radius  of  the  I-th  antenna  (meters). 

If  NJ-0,  the  rest  of  this  data  set  is  ignored. 

READ( 5 , l 8 ) ( INDJW( I ) , 1- 1 , NJ ) 

18  F0Ri*tAT(  1018) 

INDJW  -  Array  containing  the  wire  index  for  each  BOT/wire  junction. 
Each  index  must  correspond  to  a  wire  point  which  either 
starts  or  terminates  an  antenna  at  the  BOT  surface. 


READ(5, 19)(RADD(I) , I"l,NJ) 

19  FORMAT(10F8.h; 

RADD  -  Array  containing  the  radius  for  each  junction  disk  (meters). 

READ(5,20)(UXJ(I),I«1,NJ) 

20  F0RMAT( 10F8.4) 

READ(5,?)(UYJ(I),I-1,NJ) 

21  FORMAT (10F8. 4) 

READ(5,22)(UZJ(I),I«1,NJ) 

22  F0RMAT( 10F8.4) 

UXJ,UYJ,UZJ  -  Arrays  containing  the  x,  y,  and  z  components, 
respectively,  of  the  normal  vector  to  each  junction 
disk.  The  normal  vector  need  not  be  normalized  to 
unity. 

2.3.4  Aperture  Antenna  Data  Set 

This  data  set  is  read  only  by  the  BOTRA  and  BOTAC  programs  and  should  be 
deleted  when  running  the  BOTSCB  and  BOTSCM  programs. 

READ(5,23)NSA 

23  F0RMAT(I3) 

NSA  -  Number  of  slot  antennas  on  the  BOT. 

READ(5,24)(iS(K),K=l,NSA) 

24  F0RMAT( 1018) 

IS  (K)  -  Triangle  function  at  which  slot  antenna  K  is  located 
(centered). 

READ(5,25)(ZO(K) ,K“1 ,NSA) 

25  F0RMAT( 10F8-4) 

Z0(K)  -  Starting  Z  coordinate  for  antenna  K  (meters). 
READ(5,26)(Z1(K),K-1,NSA) 
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26  ?ORMAT( L0F8.4) 

Z1(K)  -  Ending  Z  coordinate  for  antenna  K  Caters). 

(ZO(K)  <  Z1(K)  for  all  K.) 

READ(5 ,27)(EO(K) ,K-1 ,NSA) 

27  FORMAT( 10F8.4) 

EO(K)  -  Constant  excitation  across  slot  antenna  K  (E0(K)  Is  complex). 

READ( 5 , 28) (TEXC(K) ,K“1 ,NSA) 

28  FORMAT( 10F8.4) 

TEXC(K)  -  Indicates  t  excitation  on  slot  antenna  K  when  TEXC(K)**1 ,0. 
Set  TEXC(K)*0.0  If  antenna  K  is  not  excited  In  the  t 
direction. 

READ(5,29)(ZEXC(K),K«1,NSA) 

29  FORMAT(10F8.4) 

ZEXC(K)  -  Indicates  z  excitation  on  slot  antenna  K  when  ZEXC(K)**1 .0. 
Set  ZEXC(K)-0.0  If  antenna  K  is  not  excited  in  the  z 
direction. 

2.3.5.  Wire/Junction  Voltage  Data  Set 

This  data  set  Is  read  only  by  the  BOTRA  program  and  should  be  deleted 
when  running  the  BOTSCB  and  BOTSCM  programs.  This  data  set  is  Ignored  If 
NW«0. 


READ(5,30)NWJV 

30  F0RMAT(I3) 

NWJV  -  Number  of  wire  and/or  junction  voltage  points. 

If  NWJV  **  0,  the  rest  of  this  data  set  is  ignored. 

READ(5,31)(IW(I),I-1,NWJV) 

31  F0RMAT( 1018) 

IW  -  Array  containing  the  NWJV  voltage  point  indices  in  the  wire 
voltage  array  to  be  fed.  The  K-th  wire  triangle  Is  fed  if 
IW(I)-K.  The  K-th  junction  is  fed  if  IW(I)«(NPW-3*NW)/2  +  K. 
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READ(5,32)(EU(  l),l-l,NWJV) 

32  F0RMAT(  10F8.4) 

EW  -  Array  containing  the  NWJV  voltages  (complex). 

2 . 3 . h  Radiation  and  Scattering  Analysis  Data  Sec 

This  data  set  is  read  by  the  BCTRA,  BOTSCM,  and  BOTSCB  programs. 

READ( 5,33) NANG , NT , P  HI I , THI 

33  FORMAT(2I3,2F8.4) 

NANG  -  Number  of  fixed  radiation  or  scattering  angles,  as  defined 
by  IPLANE  below.  NANG  radiation  or  scattering  patterns  will 
be  calculated. 

NT  -  Number  of  varied  radiation  or  scattering  angles,  as  defined  by 
IPLANE  below. 

PHI  I  -  <j>  angle  of  incident  wave  (degrees).  Used  only  in  BOTSCB 
program. 

THI  -  0  angle  of  incident  wave  (degrees).  Used  only  in  BOTSCB 
program. 

READ( 5 , 34) (ANG( I) , I-l .NANG) 

34  FORMAT ( 10F8.4) 

ANG  -  Array  of  fixed  radiation  or  scattering  angles,  as  defined  by 
IPLANE  below. 

READ( 5 , 35)  (  IPLANE( I) , 1-1 , NANG) 

35  FORMAT(IOIB) 

IPLANE  -  Array  indicating  whether  the  corresponding  element  of  array 
ANG  is  a  ip  or  6  angle.  IPLANE(I)-1  indicates  that  ANG(I) 
is  a  fixed  $  angle.  In  this  case,  <j>  is  fixed  at  ANC(I)  and 
8  varies  between  ANG1(I)  and  ANG2(I)  at  NT  equally  spac’d 
angles.  IPLANE(I)=2  indicates  that  ANG(I)  is  a  fixed  6 
angle.  In  this  case,  6  is  fixed  at  ANG(l)  and  $  varies 
between  ANGi(I)  and  ANG2(I). 

READ(5, 36)(ANG1( I ) , I“1 , NANG) 
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36  FORMAT( 10F8.4) 

ANG1  -  Array  of  starting  angles  (degrees). 

READ( 5 , 37 ) (ANG2( I) , 1-1 .NANG) 

37  FORMAT( 10F8. 4) 

ANG2  -  Array  of  ending  angles  (degrees). 

2.3.7  Near-Field  Analysis  Data  Set 

This  data  set  is  read  by  the.  BOTRA  and  BOTSCB  programs. 

READ(5,38)NTEST 

38  FORMAT (13) 

NTEST  -  Number  of  test  points  at  which  near-field  radiation  or 
scattering  is  to  be  calculated.  NTEST  may  be  set  to  zero. 

Repeat  the  following  NTEST  times: 

READ( 5 ,39) ZTEST , YTEST , XTEST 

39  FORMAT(3F8,4) 

ZTEST  -  z  coordinate  of  the  test  point  (meters). 

YTEST  -  y  coordinate  of  the  test  point  (meters). 

XTEST  -  x  coordinate  of  the  test  po^nt  (meters). 

All  of  the  above  user-specified  input  data  are  checked  for  errors  by  the 
A-STAR  codes  and  printed  on  the  line-printer  as  a  given  problem  is  executed. 

In  addition,  the  BOT  generating  curve  is  plotted  on  the  line-printer  as  a 

visual  check  of  the  BOT  input  geometry  (the  x  and  y  axes  use  the  same  scale 
for  this  plot).  The  codes  BQTAC  and  BOTRA  generate  plots  of  the  three 
principal  views  of  the  input  geometry  (these  plots  are  not  to  scale).  These 
plots  can  be  used  to  check  for  errors,  e.g.,  in  the  wire  inputs  and  near-field 
points.  Definitions  of  the  symbols  used  in  these  plots  are  given  -below: 

B  -  BOT  generating  curve  data  point 
C  -  Cap  center  data  point 
E  -  Edge  termination  data  point 
J  -  Junction  data  point 
N  -  Near-field  data  point 
S  -  Slot  antenna  termination  point 
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W  -  Wire  data  point 

+  -  Two  or  more  of  the  abcve  symbols  occupy  the  same  line-printer 

position. 

2.3.8  BOT  Inversion  Data  Sets 

The  BOT  inversion  code  (BOTINV)  is  used  when  the  'ystem  matrix  to  be 
inverted  contains  the  open  BOT  only.  The  following  input  data  are  required: 

READ( 5 , L )NM0DE , NBAND 

1  FORMAT( 213) 

NMQDE  -  Number  of  nonnegative  modes  to  be  used  for  the  inversion 
(i.e.,  there  will  be  2*NMQDE-1  total  modes). 

NBAND  -  Number  of  submatrix  diagonal  bands  to  be  used  during  the 
inversion  of  ZB0T.  HBANIW  indicates  that  only  the  diagonal 
Zm  n  submatrices  are  to  be  used  during  inversion.  NBAND«2 
indicates  that  diagonal  and  off-diagonal  Zm  n  submatrlces  are 
to  be  used,  etc.  If  NBAND  >  2*NM0DE-1 ,  the  entire  ZBq,j.  matrix 
la  inverted. 

READ(5,2)(LOAD(I),I-l  ,1IM) 

2  FORMAT( 10F8.4) 

READ( 5 , 3) (LOAD( I) , I-NM+1 , 2*NM) 

3  FORMAT ( 10F8.4) 

LOAD  -  Array  containing  the  BOT  surface  Impedance  loading  in  the  t  and 
a  directions,  at  each  triangle  peax  on  the  BOT.  L0AD(1)  through 
LOAD(Nh)  contain  the  t-directed  loadings,  and  L0AD(NM+1)  through 
L0AD(2*NM)  contain  the  z-directed  loadings,  where  NM*=(NP-3)/2 
is  the  number  of  triangle  functions  on  the  BOT.  LOAD  is 
a  complex  array. 

When  the  system  matrix  to  be  inverted  contains  the  open  BOT  with  wires 
and/or  caps  included,  the  BOTINVA  code  is  used.  The  BOTINVA  code  car.  be  used 
to  perform  the  following  types  of  system  matrix  inversion: 
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1)  Caps  can  be  added  to  a  previously  inverted  system  which  does  not 

contain  caps. 

2)  Wires  can  be  added  to  a  previously  inverted  system  which  does  not 

contain  wires. 

The  following  data  are  required  by  BOTINVA; 

READ(5,  l)NC,NPR,h'E,NW,NPW,NJ 

1  F0RMAT(6I3) 

Either  NC  or  NW  must  be  zero. 

NC  -  Number  of  end-caps  to  be  added  to  the  system. 

NPR  -  Number  of  radial  points  used  to  describe  each  cap. 

(NPR  must  be  odd). 

NE  -  Indicates  whether  or  not  the  BOT/cap  transition  term  is  to  be 
included.  NE-0  indicates  no  transition  term,  and  NE  *  0 
indicates  that  the  BOT/cap  transition  term  should  be  used. 

NW  -  Number  of  wire  antennas  to  be  added  to  the  system. 

NPW  -  Total  number  of  data  points  representing  all  of  the  wire 
antennas . 

NJ  -  Number  of  BOT/wire  junction  points. 

If  NC  *  0,  the  following  is  required: 

READ( i> , 2) (L0AD( I)  ,  1-1  ,NM) 

2  F0RMAT(10F8.4) 

READ(  5 , 3  )(I,0AD(  I )  ,  I-NM+1 , 2*NM) 

3  F0RMAT( 10F8.4) 

LOAD  -  Array  containing  the  BOT  surface  impedance  loading  in  the  t  and 
z  directions  at  each  triangle  peak  on  the  BOT.  L0AD( ] )  through 
„JAD(NM)  contain  the  t -directed  loadings,  and  L0AD(NM+1) 
through  L0AD(2*NM)  contain  the  z-directed  loadings,  where  NM“ 
(NP-3)/2  if  the  number  of  triangle  functions  on  the  BOT. 

LOAD  is  a  camplex  array. 


•VST'*1 


f<V~ 
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2.4  Parameter  Selection 

The  choice  of  most  of  the  parameters  in  the  foregoing  sections  Is 
specified  by  the  user  depending  upon  the  explicit  requirements  of  the 
problem.  For  example,  NSA  is  set  by  the  number  of  slot  antennas  on  the  BOT. 
Similarly,  BK  Is  determined  by  the  frequency  at  which  the  MM/BOT  analysis  is 
carried  out.  On  the  other  hand,  the  choice  of  some  parameters  Is  based  upon 
the  requirements  of  the  MM/BOT  theory,  as  explained  below. 

The  parameter  NMODE  is  set  by  the  length  of  the  BOT.  In  general,  the 
minimum  requirement  is  that  NMODE  >  2L/A,  where  L  is  the  (axial)  half-length 
of  the  BOT.  This  requirement  is  comparable  to  the  MM/BOR  analysis  requirement 
that  the  maximum  circumferential  modes  used  be  n  ~  flD/A,  where  D  is  the 
largest  diameter  of  the  BOR.  As  a  general  observation,  the  accuracy  of  the 
analysis  increases  and  the  spatial  resolution  of  the  surface  currents  on  the 
BOT  is  improved  as  NMODE  is  Increased.  This  trend  is  particularly  true  for 
the  edge  currents.  However,  practical  computer  main-memory  limitations 
usually  set  the  upper  limit  on  NMODE. 

The  parameters  NP,  NPR,  and  NPW  determine  the  special  resolution  with 
which  the  BOT  generating  curve,  caps,  and  wires  are  described.  In  general, 
the  maximum  separation  between  data  points  should  not  exceed  0.075X,  which  Is 
equivalent  to  0.15X  between  triangle  peaks.  As  a  general  principle,  the 
accuracy  of  the  analysis  is  improved  as  the  separation  between  data  points  is 
decreased. 

The  parameter  ZE(I)  indicates  the  z  coordinate  at  which  the  BOT/cap 
transition  triangle  function  starts  on  the  BOT.  In  general,  ZE(I)  should  be 
located  on  the  order  of  0.1 5A  from  the  end-c^n  (l.e.,  ZC(I)). 

The  parameter  RADD(I)  specifies  the  BOX/ junction  disk  radius  on  the 
BOT.  In  general,  RADD(I)  should  be  on  the  order  of  0.15X. 


Program  Dimensions 


The  A-STAR  programs  (Volume  III)  are  currently  configured  to  handle 
problems  with  the  following  set  of  parameters.  (Some  of  the  arrays  are 
overdiniensioned  in  the  listings.) 
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NMODE  <  4 
NBAND  <  2*NMOl)E-1 

NP  <  19  (17  for  a  closed  generating  curve) 

NANG  <  6 
NT  <  91 
NSA  <  20 
NC  «  2 
NPR  <  5 
NW  <  2 
NPW  <  18 
NJ  <  2 

For  a  different  sec  of  Input  parameters,  the  minimum  dimensions  required 
for  each  program  are  listed  below.  Table  1  contains  definitions  of  the 
parameters  used  as  the  subscripts  in  the  arrays  listed  in  the  dimension 
statements. 
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Table  1.  DEFINITION  OP  DIMENSION  STATEMENT  INDICES 


Parameter 


Definition 


KMODE  2*NMODE-l;  total  number  of  modes  for  the 
BOT  current  expansion 

LC  NC*NM*LR;  total  number  of  triangle  functions  on  all  end-caps 

for  one  component  of  current 

LE  NE*NM;  total  number  of  BOT/cap  transition  triangle  functions 

LR  (NPR-3)/2;  number  of  radially  directed  triangle  functions  on 

each  cap 


NMODE 


NPLOT 


NTEST 


(NPW-3*NW)/2;  number  of  triangle  functions  on  the  wires 

Scattering  and  radiation  analysis  data  set  input 

Cap  data  set  input 

Cap  data  set  input 

Wire  data  set  input 

(NP-3)/2;  number  of  triangle  functions  on  the  BOT 
BOT  data  set  input 

Number  of  points  on  the  BOT  generating  curve  (input).  If  the 

curve  is  closed,  use  (NP+2)  in  place  of  NP  in  all  definitions 

2*NP+2  *NC+NPW+2  *NSA+NTEST 

Cap  data  set  input 

Wire  data  set  input 

Aperture  antenna  data  set  input 

Scattering  and  radiation  analysis  data  set  input 

Near-field  analysis  data  set  input 

Wire  data  set  input 

Wire/ junction  voltage  data  set  input 


COMMON  Statement  Minimum  Dimensions  Contained  In  BOT  Programs 
COMMON/BOT2/NP ,BL,YB(NP) ,XB(NP) ,YBl(NP-l) ,XB1(NP-1) 
C0MM0N/B0T3/DH(NP-l),SV<NP-i),CV(NP-l) 

COMMQN/BOT5/T(4*NM) ,TP(4*NM) ,TZ(4*NM) 

COMMON/WIRE1/NPW,XW(NPW) ,YW(NPW) ,ZW(NPW) .XWl(NPW-l) ,YW1(NPW-1) ,ZW1(NPW-1) 
COMMON/WIRE2/DHW(NPW-l),DXW(NPW-l),DYW(NPW-l),DZW(NPW-l) 

C0MM0N/WIRE3/NW, INDW(NW+1) ,RADW(NW) 

C0MM0N/WIRE4/LW,TW(4*LW) ,TPW(4*LW) , INDTW(LW) 
cOMMON/JUNCl/NJ, INDJW(NJ) ,RADJ(NJ) ,RADD(NJ) 

COMMON/ JUNC2/TJ(2*NJ),TPJ(2*NJ),INDTJ<NJ) 

COMMON/ JUNC3/XJ(NJ)  ,YJ(NJ)  ,ZJ(N.J) 

COMMON/JUNC4/UXJCNJ) ,UYJ(NJ) ,UZJ(NJ) 

COMMON/ JUNC 5/ UXJ 1 ( NJ ), UYJ 1 ( NJ ), UZJ 1 ( NJ ) 

COMMON/ JUNC6/UXJ2(NJ) ,UYJ2(NJ) ,UZJ2(NJ) 

COMMON/SLOTl/NSA,IS(NSA) ,ZO(NSA) ,Z1(NSA) 

COMMON/SLOT2/EO(NSA) .TEXC(NSA) .ZSXC(NSA) 

COMMON/  CAP1/NC ,  XC ,  YC ,  ZC(  NC) 

COMMON/CAP2/NPR , RHOC ( NPR) , RHOC l(NPR-l) , DRHOC ( NPR-1) 

COMMON/ CAP3/TCR(4*I,R) ,TCT(4*LR) ,TPCR(4*LR) ,TPCT(4*LR) 
COMMON/CAP4/RC(NP),RCl<NP-l),AC(NP-l),CPC(NP-l),SPC(NP-l) 

COMMON/EDGI /NE,ZE(NE) ,ZBE(2*NE) 

COMMON/EDG2/TCE(2*NE) ,TPCE(2*NE) ,TBE(2*NE> ,TPBE(2*NE) 

COMMON '  PLOT  l  /  NPLOT ,  XPLO'i'(  NPLOT  )  ,  YPLOT(  NPLOT)  , ZPLOT( NPLOT) ,ISYM( NPLOT) 


BOTZSN  '.mum  Dimensions 

COMPLY  Z(LS*LS)  ,G((NP-l)*NP/2)  ,GB((NP-l)*NP/2) 

DJMFKt:  '.OS  TWGHT(NM)  ,ZWGHT(NM) 

BGTZSW  Ein Lnv.ra  Dimensions 

JMI’Lt «  Z(I,S*( LW+N J )  ) , G 1  (  ( NP- 1 ) * ( NPW- 1 )  )  , G 2 (  ( NP-1  )* ( NPW- 1 ) ) 

BOTZWW  X. iLraiUt  Dimension 


COMPLEX  Z(n,W+NJ)**2) 
COMMON/ WI RE 1 / NPW , X( NPW , 3 ) 
COMMON/ A7ri;2/WL(NPW-1) 
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C0MM0N/WIRE3/NW, INDW( NW+1 ) , RADW(NW) 

C0MH0N/WIRE4/TWIRE ( LW , A ) , TPW( LW , 4 ) , UW( NPW- 1 , 3 ) . 
C0MM0H/JUNC1/NJ ,RADJ(NJ) ,RADD(NJ) 

COMMON/ JUNC2/TJUNC(NJ, 2) ,TPJ(NJ,2) 

COMMON/ JUNC3/XJ( N J , 3 , 3 ) , WLJ ( NJ , 2 ) 

COMMON/ JUNC4/UJ( NJ , 2 , 3) 

COMMON/ JUNC5/URT(  N.T ,  3  ) 

COMMON/ JUNC6/UFZ{NJ, 3) 

DIMENS  ION  INDJW(  NJ  ) ,  INDT  J(  NJ  ) , UXJ ( NJ ) , UYJ( NJ ) , UZ J( N J ) 
B0T2SC  Minimum  Dimensions 

COMPLEX  Z(LS*(2*LC+LE)),G2((NP-l)*NC*(NP-l)*(NPR-l)) 
COMPLEX  GIE( (NP-1 )*NE*(NP-1 )*2) ,G2E<  <NP-i*NE*(NP-l )*2) 

BOTZCC  Minimum  Dimensions 

COMPLEX  Z(4*LC*LC),GS((NP-1)*(NFR-I)> 

DIMENSION  TWGHT(NM*LR) ,RWGHT(NM*LR) >EWGHT(NM*NE/NC) 

BOTZCW  Minimum  Dimensions 

COMPLEX  Z((2*l,C+LE)*(LW+NJ)) 

BOTINV  Minimum  Dimensions 

COMPLEX  Z(K1),ZI(K2),  LOAD(K3) 

DIMENSION  WGHT(K3) 

DIMENSION  NZ(K4) 

COMMON  \M,  JK(4) ,LR(K5) 

where  KI  through  K5  depend  on  NBAND  and  ars  defined  below. 


NBAND 

1 

<2*NM0DE-i 

>  2*NMODE-l 

Kl 

0 

LS2*  { (2*NM0DE-1) 

* ( 2*NBAND- I ) - ( NBAND- i )*NBAND } 

|lS*(2*NMODE-1)} 

K2 

LS2 

LS2*(2*NMODE-l) 

LS2 

K3 

LS 

LS 

LS 

K4 

0 

2*NM0DE-I 

0 

K5 

LS 

LS 
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LS* ( 2*NM0DE-1 ) 

BOTINVA  Minimum  Dimensions 


The  program  dimensions  depend  on  the  type  of  matrix  Inversion  being 
performed.  The  four  possible  types  along  with  the  corresponding  dimensions 
follow: 

1)  Addition  of  wires  to  an  open  BOT 

COMPLEX  PI( (LS*KMODE)2) ,Q(LS*KMODE*(LW+NJ) ) ,R(LS*KMOD£*(LW+NJ)) ,S( (LW+NJ)2) , 
YI(max(LS*LS, (LW+NJ)2)) 

COMPLEX  Wl(max(LS*KMODE, LW+NJ) ) ,W2(LW+NJ) 

COMPLEX  L0AD(1) 

DIMENSION  WGHT(l) 

2)  Addition  of  caps  to  an  open  BOT 

COMPLEX  PI( ( LS*KMODE)2) ,Q(LS*KMODE*( 2*LC+LE) ) , R(LS*KM0DE*(2*LC+LE) ) , 
S((2*LC+LE)2),  YI(max(LS*LS,(2*LC+LE)2) 

COMPLEX  Wl(max(LS*KM0DE,2*LC+LE)) ,W2(2*LC+LE) 

COMPLEX  L0AD(2*LC+LE) 

DIMENSION  WGHT ( 2 *LC+LE ) 

3)  Addition  of  caps  to  an  open  BOT  with  wires 

COMPLEX  PI ( ( LS*KMODE+LW+N J ) 2 ) , Q ( ( LS*KMODE+LW+NJ ) * ( 2 *LC+LE ) ) , 

R( ( LS*KMODE+LW+N J ) * ( 2*LC+LE )  )  ,  S  (  ( 2*LC+LE) 2 ) , 

YI (max( LS*LS , ( LW+N J ) 2 , ( 2*LC+LE) 2 ) ) 

COMPLEX  W 1 (max ( LS*KMODE+LW+N J , 2 * LC+LE ) ) , W2 ( 2  * LC+LE ) 

COMPLEX  LOAD ( 2*LC+LE ) 

DIMENSION  WGHT(2*LC+LE) 

4)  Addition  of  wires  to  a  BOT  with  caps 

COMPLEX  PI((LS*KMODE+2*LC+LE)2)  ,Q( (LS*KM0DE+2*LC+LE)*(LW+NJ) )  , 

R( ( LS*KM0DE+2*LC+LE ) * ( LW+N J ) ) , S ( ( LW+N J ) 2 ) , 

YI (max( LS*LS , ( LW+N J ) 2 , ( 2*LC+LE ) 2 ) 

COMPLEX  W1 (max( LS*KM0DE+2*LC+LE , LW+NJ ) ) , W2 ( LW+NJ ) 

COMPLEX  L0AD(1) 

DIMENSION  WGHT (1) 

BOTAC  Minimum  Dimensions 

COMPLEX  Y(MAX( LS*LS , ( LW+NJ ) 2 , ( 2*LC+LE ) 2 ) ) 

COMPLEX  TAB(NSA*NSA) ,ZAB(NSA*NSA) ,WAB(NJ*NJ) 
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BOTRA  Mi  i: mum  Dimensions 


% 

I 

*> 
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COMPLEX  v:if  LS)  ,VW(LW+NJ) 

COMPLEX  GT(is.'"',GP(NT) 

COMPLEX  C8( LS*XMODE) ,CW(LW+NJ) ,CC(2*LC+LE) 

COMPLEX  Y(max(  LS*LG ,  ( LW+N J  )  2 ,  (  2*LC+LE )  2 > 

COMPLEX  RBT(LS) ,RBP(LS>  -RWT<  LW+NJ) , RWP(LW+NJ) , RCT( 2*LC+LE) ,RCP(2*LC+LE) 
DIMENSION  THR(NT) , PHIR(NT) 

DIMENSION  ANG( NANG ) , I P  LANE ( NANG ) , ANG 1 ( NANG ) , ANG2 ( NANG ) 

SUBROUTINE  VWIRE 
DIMENSION  IW(NWJV) 

COM°LEX  EW(NWJV) 

SUBROUTINE  NEARB 

COMPLEX  GT(NP-1),GZ(NP-1),G1T(NP-1),G1Z(NP-1),H1T(NP-1) 


BOTSCM  Minimum  Dimensions 

COMPLEX  STT(NT) ,SPP(NT) ,STP(NT) ,SPT(NT) 

COMPLEX  CBT(LS)  ,CBP(LS)  ,CWT(  LW+NJ  ) ,  CWP  (  LW+N  J  )  ,  CCT  C  2*L.C+LE  ) ,  CCP  (  2*LC+I.E  ) 
COMPLEX  Y (max( LS*LS , ( LW+NJ ) 2 , ( 2*LC+LE) 2 ) 

COMPLEX  RBT( NT*LS ) , RBP( NT*LS ) , RWT ( NT*( LW+NJ ) ) , RWP ( NT* ( LW+NJ ) ) , 

COMPLEX  RCT( NT*( 2*LC+LE ) ) , RCP( NT* ( 2*LC+LE) 

DIMENSION  THS(NT) ,PHIS(NT) 

DIMENSION  ANG(NANG) , IPLANE(NANG) ,ANG1 (NANG) , ANG2(NANG) 


BOTSCB  Minimum  Dimensions 

COMPLEX  STT(NT) ,SPP(NT) ,STP(NT) ,SPT(NT) 

COMPLEX  CBT(LS*KMODE) ,CBP(LS*KMODE) ,CWT( LW+NJ) ,CWP(LW+NJ) ,CCT(2*LC+LE) , 
CCP(2*LC+LE) 

COMPLEX  Y( max ( LS*L  S , ( LW+NJ ) 2 , ( 2*LC+LE ) 2 ) 

COMPLEX  RBT(LS) ,RBP(LS) ,RWT( LW+NJ ) , RWP ( LW+NJ ) , RCT ( 2 * LC+LE ) f RCP ( 2*LC+LE ) 
DIMENSION  THS(NT) , PHIS (NT) 

D IMENS ION  ANG ( NANG ) , IPLANE ( NANG ) , ANG i ( NANG ) , ANG2 ( NANG ) 

SUBROUTINE  NEARB 

COMPLEX  GT(NP-l) ,GZ(NP-1) ,G1T(NP-1) .GlZ(NP-l) ,H1T(NP-1) 
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2.6  Sample  Problems 

In  this  section,  four  sample  problems  are  considered  to  Illustrate  the 
use  of  the  MM/BOT  algorithm.  Sampxe  Problem  i  exercises  all  main-line 
programs  given  in  Figure  1  for  BOTe  having  aperture  (slot)  antennas.  The 
coordinatr  generation  and  specification  of  antenna  location  and  feed  are  i 

demonstrated.  The  inputs,  outputs,  and  selected  intermediate  results  are  j 

reproduced  here  in  detail  to  provide  a  check  case  for  the  proper  functioning 
of  the  codes.  Problem  2  demonstrates  the  radiation  analysis  for  a  monopole 

t 

antenna  attached  to  a  BOT.  Inclusion  of  the  junction  effects  exercises  the  * 

junction-related  parts  of  the  impedance  generating  routines  together  with 

those  for  the  wire  representation  of  the  monopole.  The  second  and  third  parts 

of  Problem  2  demonstrate  the  analysis  for  two  monopoles,  one  active  and  the 

other  passive,  and  a  loop  antenna.  In  Problem  3,  the  radiation,  near-field, 

and  coupling  analyses  are  carried  out  for  an  active  and  passive  monopole 

mounted  on  the  trailing  end  of  an  asymmetric  wing  section.  Problem  A  utilizes  j 

the  scattering  analysis  routines  for  a  closed  cylindrical  body.  In  this 

example,  the  procedures  are  demonstrated  for  inclusion  of  the  edge  transition 

region  between  the  caps  and  the  BOT  surface. 


2.6.1  Problem  1:  Aperture  Antenna  on  BOT 

Consider  a  right-circular  cylinder  of  2.761  length  and  0.2161  radius  with 
an  embedded  ^-polarized  aperture  antenna  at  $  »  90°.  The  aperture  is  fed 
uniformly,  subtends  a  45°  opening,  and  is  2.07a  long  in  the  axial  direction 
(see  Figure  5a).  For  simplicity,  consider  that  the  BOT  is  uncapped.  (The 
radiation  pattern  for  the  capped  body  is  substantially  the  same  as  for  the 
present  case.) 

a)  Calculate  the  power  gain  patterns  in  the  horizontal  (<{,  =  0,  l-“0°)  plane 
and  the  roll  (9  ■  ±90°)  planes  in  the  6  and  $  polarizations. 

b)  Compute  the  currents  on  the  cylinder  surface. 

Solution  -  The  calculations  were  carried  out  at  10  MHz  (1  »  00  ra),  using  four 
modes  (NM0DE=4)  to  represent  the  BOT  currents.  The  cylinoer  was  represented 
by  N?»17  points  around  the  circumference. 

B0TSEG  was  executed  in  a  time-share  mode  (Figure  5b)  in  order  to  obtain 
the  BOT  generating  carve.  Figure  5c  lists  the  input  data  used  to  execute  the 
programs  BOTZSS  and  BOTRA.  For  reference,  the  variables  of  the  data  set  are 
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labeled.  The  aperture  coincided  with  the  second  triangle  function  (IS  *  2), 
and  only  one  triangle  function  was  used  to  span  the  aperture  (NSA  ■  1).  If  a 
nonuniform  aperture  excitation  is  desired,  then  more  triangle  functions  should 
be  used  to  span  the  aperture,  each  with  a  different  EQ.  BOTINV  was  executed 
using  NMODE-4,  NBAND-14,  and  zero  surface  impedance  loading. 

Partial  outputs  from  BOTZSS,  BOTINV,  and  BOTRA  are  shown  in  Figures  6,  7, 
and  8,  respectively.  The  radiation  power  gain  for  the  slotted  cylinder  for 
the  vertical  plane  normalized  to  an  isotropic  radiator  is  summarized  in  Figure 
8a.  (The  comparison  of  these  results  with  the  MM/BOR  analysis  is  given  in 
Figure  8  of  Volume  I.)  Partial  output  of  the  currents  on  the  cylinder  is 
plotted  in  Figure  8b. 

The  foregoing  calculations  were  carried  out  at  10  MHz.  If  the  dimensions 
of  the  body  (BOT)  and  the  antenna  are  given  Initially  in  terms  of  wavelength, 
any  convenient  frequency  can  be  chosen  in  the  setup  procedure  for  carrying  out 
the  computations.  If  the  data  are  to  be  compared  with  range  measurements  at  a 
given  frequency,  for  ease  of  data  Interpretation,  the  calculations  also  are 
done  at  that  frequency. 


Figure  5a.  Slotted  cylinder  for  Problem  1. 
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Figure  5b.  Execution  of  BOTSEG  for  Problem  1 ,  with  triengle  functioni  numbered  on  the  plot. 
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Figure  6a.  Partial  output  of  BOTZSfi  for  Problem  1 . 


3u 


2.4798 

2.466' 


1  Hu  \ 

NJMBf  R^OP^nUQHATIQNS^JND  MCRAGt  NUMBER  OP  POINTS  PER  INTEGRA  f  ION  ISEl.7  A4U  NUN»*E«.F  URNS.  RESPECT  tvELT) 
l  Sill  •  Stl) 

,1929E*C5  -.45271*05  ilSl4C*0S  .3U6if.ps  TiOp.  .16181*05  5445,  .156TC-0S  »1M.  J709E* 

•  'fill*!?  j|wtt»»s -"•»}{  !*j{  •  it tick.  *♦,*.  .ism* 

Imi  :  U$|:!i  -fill  :f|  li!lt£*°‘5  »*•••  •>»'«• 

6190.  ,  'OOf*Of  5449.  .!>•»  Of  TIoo,  .  616E*05  .15141*05  .10611*05  .  192H*05  -.*5271* 

,lsI*E*05  .V  611*05  7100.  .Ul»;*ps  544*.  .  5*71*05 

*}; :jl*  ::  j  TttS:  : tt«:  f  nos-  •,4“t-°4  •14“t'44  -“M£* 

Mop,  ,uUF»of  liki.  .  1  So t  •3?  oiSi.  •  T65c*5?  4*4*.  .i66Ji«05  moo. 


14£*05  .ICGI  «OS  . 

h'-04  ! 


,1S|4£*0S 

5*i«. 

.19291*05 


l  S Mi  -  S«T» 

jioOJE-OP  'jlJ5°C'09  ^ 

#:|i  :|||-"  i’iijicii 

•fflt:  -hi!:  •«!( 

I  Mr I at.., 

.UT41-09  -t2134£-Q9  .1)491*05 

:llil;  -lift::; 

‘.SUu-W  ,62QTE-I0  *11491*05 

:M:«  L.„  :i  !■••• 

*3359.  Uao.  *41*n*os 

12Q.9  -2695.  •i|04C-09 

Mi  m...  §\m 

Jh:i  lib 
iF'liE  m 
iriilir; 

•  ins  !  $  *«•«» 

:  i  Ml  'ill  •  .%!»•.» 

;§tw*oI  -.inif'di  litof 


-»n«.  im, 

iili.  -JJ5J, 

.?fUk-,.  10 

1»|;  :i  M’u 

iii 
lip."'  ’il 

4ll4.  * 13971-09 

-4ll*.  3359. 

-Iii;,.,.  ft#- 

-siil; 


-•*40. 

-Sit:: 

•iill: 


Iii*-  : 


■jifSii-n 


:»Uc-io 

..fiiJt-o-J 


IS  P 

:•  i!i!:i  i!  ‘  " 

IM,  ,«rm«  os 

:•!!;:!!  Hi*-* 

im.  ,miH'05 


m 


544*. 

•1S67C*0S 

6190. 

■  1  7091*05 

7100, 

.U1BC-0S 

S449. 

,1567E*05 

• 15161*05 

• 3061E • OS 

7100, 

•  1 61  BE  *  05 

.  1«29£*U5 

•,*S2TC«0S 

. lS14t*0S 

•  306 1C  *05 

•15141*05 

. J061E  *  OS 

.  192H*0S 

-,•5271*05 

7100. 

•1910E-0S 

.  IS14E.0S 

.  1061E  *05 

SM«, 

•1SP7E«0S 

7100. 

. 16101*05 

6l9d. 

.1709E-0S 

S44*, 

.15671*05 

121.0 

•269f'. 

-.1940E-11 

,06191-10 

3359. 

—4040 . 

320.9 

-2695. 

•  lUOl.Qb 

•4314. 

31S9. 

-4B40. 

•1J9TE-09 

-.2129E-09 

•1I«9E*0S 

•4314. 

••ll*9£*05 

•  314. 

•  U64E-C9 

-,2599£-0? 

-3JS9. 

40<  1, 

••U*H*OS 

4314. 

-121.0 

29  9S. 

•3119. 

40*0  . 

-.UOOfOV 

-.U1SE-10 

-320.9 

2*95. 

120.9 

-299i. 

•9S04E-10 

-,10941-10 

33S9. 

-» 040. 

321.0 

-2695, 

•11 49E-0S 

-43.4, 

3359. 

-4040. 

UG4E-09 

.3lklE-lO 

,114*1405 

•4  ll  4  . 

-  •  1 1 49£ • OS 

♦41* 

• •4903C- 1 0 

•09221*10 

on*. 

4040. 

-  •  | ;49t  *05 

4114. 

•  120,9 

249%. 

- J3S9, 

40*0  . 

•JOOSt-IV 

••4jfr4i- | | 

-321.0 

•  t*.  5, 

19. Td. 

-.1904E*OS 

•4023, 

-•10601*05 

,2?2tE*0‘J 

-• 3. 1 01 • OS 

1918. 

-.19041*05 

•srioi.o^ 

•  70S, 

•  2720k  *05 

-.i*m*os 

•  M9JC.0S 

•SJ14E-0? 

•5?jOC*G5 

6706. 

.SM9C*ot» 

670S. 

,T19ii«0S 

,51141*05 

.iU  8E*0S 

-«  32  1  BE  *0S 

•  57 1 9E  *05 

6706. 

1910, 

-. 3004(«OS 

.27201*05 

-*121il*06 

-402/. 

•• 3BSOE«OS 

3938. 

-.39041*9* 

Fiyure  6*.  Concludfd. 


urtUUC  NtjlNO 

i*  *  i* 

I  And  l  6UN*4CL  tMr*EOANCl  wOAOINO  tCU4Ml&X! 

□ . Qutlg 


NMODE  and  NBAND  ate  user  specified 


o.nono  O.uoQu  *)  •  0 0 0  0  Q.QCQQ  U.OuUu  u.nDUu  U.OQGO  0*0000  O.QUOU  O.uUOO^ 

0*0000  O.OOOU  0.4  »>00  C.OOOd  O.uPUO  0*0000  I 

r  Surface  loadings  are  user  specified 

0.0000  0.0000  0 « QUOU  Q.UUOuJ 


0 .  <1000  O.UOOO  0.0  000  0.0000  o.ooov 

0 . 0 uUU  0.0000  - -  - -  “ 


0.0000  0.0000  Q.OUOu 


OOg 


I  «L  MINIMUM  dOllNV  A«WAY  DIMENSIONS  NH  1*1$  M«U8t.EM  A«t  AS  7UwLO«S» 
l  L 1  tO*C>i«&NT  NZ  t> 

12*>4a  256  16  0  IW 

fO UL  MATRIX  iMVtHSJWN 


••301 SE -OS 
•.2*771-06 
-•20a?«~05 
.*A6?e«06 
-.26421-0$ 
-•22771-06 
-.JQ14E-05 


-•10771-05  .27331-05 

-•74961-06  -• 30 l 5f -05 
-•70911-00  -•22771-06 
•lOjTC-05  -.211621-05 
•1*441-05  .4*671-04 

•  10J/1-05  -•*0421-01 
-./090t-00  -.22711-06 


/  sstt 
-3,-3 


.27J31-06 

-•74961-06 

-.30151-05 

.ACS2E-06 

•37951-10 

-.43371-19 

.61661-06 

,48541-05 

-•40521-04 

•132/1-06 

••37461-09 

-.61001-04 

.•a  JOE- 1ft 

•14121-04 

-.13271-04 

vSSZt 

•13271-06 

•11661 *|7 

.23721-1 7 

-3.-3 

•61IPC-06 

-•16121-04 

.13271-06 

•4U52C-06 

•57461-65 

•61001-05 

•22071-17 

••605*1-05 

•40521-06 

1 

.40521-06 

-.30691-17 

•29021-11 

j 

•61001-06 

.605*1-05 

-.40621-06 

.132/1-06 

-•97461-05 

-. 61001-06 

I 

.16101-17 

.14121-04 

-.13271-06 

vsttz 

•  13271-06 

-•30221-17 

.60901-11 

y-3,-3 

.61000?* 

-,14121-04 

< 1 3271-06 

•6QS21-06 

•97-61-05 

•61011-06 

'.16671-17 

-•40941-05 

.40521-06 

•17201-06 

•11641-05 

-, 14201-04 

-.31571-06 

•77301-06 

•17211-04 

•61101-06 

.25111-05 

-.31571-04 

••63911-06 

•26631-05 

•01101-04 

V«ZZ 

.61101-06 

.407 | t-05 

-•63011*04 

y-3,-3 

.3167004 

•26631-05 

•61101-04 

•17201-06 

•25111-06 

-•31571-04 

• 

••16201-04 

•77J01-06 

•172OE-04 

1 

Figure  7.  Partial  output  of  BOTINV  for  Probtam  1  (admittance  matrices). 
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Fiji  ra  8b.  Output  of  BOT  currant*  for  Problam  1. 


i  itii  .j'-  •„  r.  *»,-  * 


SI*  *-  ■’  'v  *:  "  ■  < 

•  •  f  ^  r  .  *  'V-*1  i j  fc  .  'uaa\  tt  . 

v. * :  *  ** 

f-‘-  /  ’  -'T-r.v.'  "  ■ 


t..  «2,-:  - . 


sB-*.4- 


* 

V 


i 


ij 

'i 


* 


I 


I-DIREC1E0  CURRENTS  ON  TRIaNGUE  FUNCTION  1 


Z-OlRECTEU  CORRtNTS  ON  TRIANGLE  FJNCTTUN 
*  HAGN I  T Ut»E 


*  PHASE 

IS!  , 

126.0 

'm 

!i:§ 
0.0 
-18.0 
-36,0 
-5A.0 
-72.0 

•1.  2.0 

■IfO.O 


Axial  distribution 
of  currents  on  80T 
along  strip  sub¬ 
tracted  by  triangle 
function  1. 


Figure  8b.  Concluded. 


37 


■WU"" 


0.9009  1.2960  2. 5920  3.9890  $.19*0  6.6809  T.7760 

36  POINTS  PLWTIEO 


Figure  8c.  Output  of  oaomatry  for  Problem  1. 
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2.6.2  Problem  2:  Wire  Antenna  on  BOT 

Consider  a  right-circular  cylinder  of  1.01  length  and  0.0751  radius,  vlth 
a  quarter-wave  monopole  mounted  at  the  cylinder  mid-point.  The  antenna  is 
base-fed. 

a)  Calculate  the  power  gain  patterns  in  the  horizontal  (41  «  0,  180°)  plane 
and  the  roll  (6  »  *  90°)  planes  in  the  8  and  $  polarization.  Use  a  capped  BOT 
representation  for  the  cylinder,  including  a  junction  region,  but  exclude  the 
edge  region  between  the  caps  and  the  BOT. 

b)  Repeat  the  foregoing  calculations  with  a  parasitic  quarter-wave  mono- 
pole  added,  mid-point  between  the  active  element  and  the  end  of  the  cylinder. 

c)  Repeat  the  calculations  in  a)  replacing  the  monopole  with  a  loop,  fed 
at  the  midpoint  of  the  cylinder,  with  an  off-set  of  0.037X  and  a  total  length 
of  0.25X. 

Solution  -  The  calculations  were  carried  out  with  X  »  0-508  m,  uriug  four 
modes  (NMODE-4)  to  represent  the  BOT  currents.  The  cylinder  was  represented 
by  NF  »  17  points  around  the  circumference.  Two  end-caps  were  included  (NC  «* 
2)  with  five  points  in  the  radial  direction  (NPR  -  5).  Figures  9a-9c  list  the 
input  data,  Including  wire  coordinates,  for  cases  2a-2c,  respectively. 

The  three  problems  were  solved  by  first  executing  the  programs  B0TZSS, 
BOTZSC,  and  BOTZCC  using  the  data  file  In  Figure  9a  (see  Figures  10-11  at-  13- 
15  for  partial  outputs).  The  open  BOT  system  matrix  was  generated  by 
executing  BOTINV  with  NMODE  “4.  NBAND  "  14,  and  zero  surface  loading  (see 
FLgure  12  for  a  partial  output).  The  closed  BOT  system  matrix  was  generaced 
by  executing  BOTINV A  with  NC  -  2,  NPR  -  5,  NE  »  0,  NW  -  0,  NFW  ~  0,  NJ  -  0, 
and  zero  surface  loading  (see  Figure  16  for  a  partial  output).  Next,  the  wire 
impedance  matrices  for  cases  2a-2c  were  generated  by  executing  BOTZSW,  BOTZCW, 
and  B0TZWW  using  the  data  ‘"lies  listed  in  Figures  9a-9c,  respectively  (see 
Figures  17-20,  23-25,  and  28-30  for  partial  output).  The  system  matrix  for 
the  closed  BOT  with  wires  is  obtained  by  executing  B0TINVA  for  each  case.  The 
parameters  for  case  2a  are  NC  ■  0,  NPR  “0,  NE  »  0,  NW”  1,  NPW  -  9,  and  NJ  » 

1  (see  Figure  21  for  a  partial  output).  The  parameters  for  case.  2b  are  NC  - 
0,  NPR  -  0,  NE  -  0,  NW  -  2,  NPW  -  18,  and  NJ  -  2  (see  Figure  26  for  a  partial 
output).  The  parameters  for  case  2c  are  NC  ”  0,  NPR  “0,  NE  “  0,  NW  -  1,  NPW 
■*  13,  NJ  ”  2  (see  Figure  31  for  partial  output).  Once  the  system  matrix  for 


the  closed  BOT  with  wires  is  obtained,  the  final  results  are  calculated  by 
executing  BOTRA  with  the  appropriate  data  fij.e  from  Figures  9a-9c  (see  Figures 
22,  27,  and  32  for  partial  output).  In  each  case,  the  first  junction  point  is 
fed.  The  index  of  this  point  in  the  wire/ junction  voltage  array  is  (NPW- 
3*NW)/2+l. 


.12294118*02  BK 
4  20  14  NMODE.  NPT  NBAND 

17  NP 

•  0000  .014b  .0269  • 0  352 

-.0269  - • 0352  -.0381  -.03S2 

.0000  .0029  .0112  .0235 

.0650  .052,*  .0381  .0235 

.2540  BL 

2  5  0  NC.NPR.NE 

.0381  .0000  XC.YC 

.2540  -.2540  ZC 

.0000  .2500  .5000 


.0361 

-.0269 

.0301 

.0112 


7500  1.0000 


.0352  .0269  .0146 

•.0146  .0000 

.0527  .0650  .0733 

.0029  .0000 


.0000 

.0762 


Cap  data  sat 


RHOC 


■ . U 1 46  YB' 
.0733  XBl 


NW.  NPW,  NJ 


►  BOT  data  set 


.0381 

•  0381  .0381 

.0381 

.0361  .0381 

.0381 

.0381 

.0381 

xw 

.0381 

.0540  .0698 

.0857 

.1016  .1175 

.1334 

.1492 

.1651 

YW 

.0000 

.0000  .0000 

.0000 

.0000  .0000 

.0000 

.0000 

.0000 

ZW 

1 

INDW 

.0011 

RAOW 

1 

INDJW 

.0450 

RADD 

•  oooo 

UXJ 

1.0000 

UYJ 

.0000 

UZJ 

J 

0 

1 

4 

NSA 

1 

i 

>  Wire/junction  voltage  data  set 

1 .0000 
4  37 

.  OOOu 

i 

l 

90.0 

0.0  -90.0 

90.0  j 

1 

0,0 

1  1 

0.0  0.0 

2  1 
-90.0  i 

>  Radiation  angle 

180.0 

180,0  180.0 

90.0 

\  Wire  data  set 


NTEST 


Figure  9s.  Input  data  for  execution  of  Problem  2a. 
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Figure  Ob.  Input  data  for  execution  of  Problem  2b. 
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Figure  9c.  Input  dftti  for  execution  of  Problem  2c. 
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Figure  10.  Output  of  EOT  input  dau  for  Problem*  2»  2c 
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Figura  11.  Partial  output  of  BOTZSS  for  ProUamt  2a-2c  (Z,m). 
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Figure  12.  Partia1  output  of  BOTINV  for  Problems  2a  2c 
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Fighr*  13.  Output  of  tap  input  d«t t  lot  Probl.mt  2,-2c. 


NJMtffR  0 

2592 

l  S<f»  - 


EG"M  IU*S  A  MO  AVERAGE  NUUJc.*  OF  POINTS  OEM  INTEGRATION 

0.0 


-.fom-ot: 
*.12701-01 
-.41SSE-02 
.♦OOTg-oi 
-.4'55E*02 
*.l270E-0l 
-.  70111*02 
-.51611*02 

.  «zm- 02 
.12821*01 
.J667E-02 
-•%090E*0l 
.36671-02 

•12021*01 

•72271-02 

.64111*02 

. J904E-g* 
.15191*01 
.43041*01 

•29121*14 

-.*3d*E-©l 

-•i6in-oi 

-• 39841-02 
-•19531-14 
■56y41*02 
•1?*2C-01 
.45-31-01 
.80*31*14 
-.45481-01 
-.17421-01 
-.66041*02 
.38171-14 


».6Z2bL-gz 
-.17611-01 
-.?47S:-gj 
.80931*01 
-.24751-01 
-.17611*01 
-.622*1-02 
-. *3041-02 
-.62B5E-U2 

•.it^te-oi 

27201*01 

.85991*01 

-.27201-01 

-.1757E-G1 

-.42851-02 

-.15741-02 

-.87851*02 
*■20221*01 
3024**01 
*•81671*14 
•302*1*01 
•20221*0 1 
.87051-02 
-•2762E*)4 
•87701-02 
•27321-01 
v  30 1 91-01 
•72501*14 
-.3010**01 
•.20321*01 
-.a770£-02 
.2*651-14 


~.Slbftt-02 
-.70051-02 
-.12731-01 
*• 4 lb2E*02 
.4io3E-o; 
-.4;*2L*ft2 
-.12731-01 
-.7ooSt-02 

•64101-02 

.72211-02 

.128*1-01 

•36751-02 

-.40951-01 

.36751-02 

•12041-01 

.72211-02 

- .15681*14 
•39841*02 
•15221-01 
•43791*01 
•19411*14 
*•43791*01 
-•15221*01 
*.1984£*02 
•3471E*14 
•56001*02 
>  1  7461*01 
•45431.01 
•70121*14 
“*4  j43fc**Qi 

*•17451*01 

-.56001-02 


• .4)091-02 
-.62241*02 
*. 1 764E-0 1 
-.24751*01 
.89091*01 
-.24751-01 
-.17641-01 
-.62241-02 
•• I 5751-02 
-.42001-02 
-.17601-01 
-.272IE-01 
.86061-01 
-.2?2lt-0i 
-.17601-01 
-.42801-02 

-.32321-14 

-.87761*02 

-.20251-01 

-•30191-01 

-.62241-14 

•30191-01 
•20251*01 
.87761*02 
.22191*14 
•87691*02 
.20351-01 
•30151-01 
•5)001* ) 4 
••30151*01 
-•20351*01 
*,87691*02 


riourt  14.  Partial  output  of  BOTZSC  foi  Probtamt  2a-2c. 
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Figure  15.  P«rtUI  output  of  BOTZCC  for  Probl.mt  7a-2c. 
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Figur*  16.  Concluded. 
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Figure  16.  Partial  output  of  BQTINVA  for  Problems  2a-2c  when  caps  are  added  to  the  system  matrix. 


WIRE  C30R0!n*TE* 

I.  X* 


xe 

re 

l* 

0361 

•  0361 

•  Ub4Q 

o.oooo 

oaai 

0.0000 

0361 

.0696 

o.oooo 

in! 

.U6S7 
•  1 616 

0.0000 

c.0000 

036, 

•  1175 

o.oooo 

0361 

•  1334 

0.0000 

0301 

•  1*92 

o.tooo 

0  361 

.1651 

0.0000 

JUNCTION  PXRXHEIERS 
1 J  RXOO  OXJ 


uxoo  jij  UTJ  Ulj 
.0*50  U.OOOU  1.0000  0.0000 


Figure  17«.  Output  of  wire  input  dot*  for  Problem  2a. 
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Figure  17b.  Output  of  wire  input  dot*  for  Problim  2b. 
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Figure  17c.  Output  of  wir#  input  deu  for  Problem  2c. 
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Figure  13.  Partial  output  of  BOTZWift  for  problem  2a. 
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Figure  19.  Ptrtitl  output  of  30TZSW  for  Problem  2» 
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Figure  20.  Partial  output  of  BOTZCW  for  .roblem  2a. 
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Figura  24.  Partial  output  of  BOTZSW  for  Probiei.i  26. 
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Figure  25.  Partial  output  of  BOTZCW  for  Problam  2b. 
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Figure  27.  Pertlil  output  of  BOTRA  for  Problem  2b. 
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Figure  28.  Partial  output  of  BOTZSW  for  Problem  2c. 


l  cm  - 

-.44166 

-.14126 

-.2100E 

•88236 

-.49036 

-.4896E 

I5U7E 

:;t  !ttf 

-.34526 


*  I  HE 

-0?  -.77*36-02  - 

-14  .75j2E-15  - 

-03  .44646*02  - 

-04  -.62-:3E'02 

:ii  =:|S3?i=Si  : 

-0?  -• 14016-02  - 
-02  .17266-02 

-02  -.99906-02  - 
-14  .2311E-15  - 

*02  -.72946-04  - 


.47316-02  -.4J22F-04 
-.468lf-02  -.10976-01 
-.40176-14  -.30406-14  - 


-.27,876- 

*• b|42E- 
.72546- 


02  .83036-03 
02  -.12496-02 
0?  -.64036-02 
14  -.19466-14 
02  .10276-02 
02  -.22916-02 


..S649E- IS 
-.44/46-02 
-.88236-04 

-.56576-02 

-.*7316-02 

.20066-14 

■.*4756-02 

-.3668E-02 

-.*0176-14 

-.91166-02 

-.72546-02 


.10046-14 
-. >0306-02 
.02236-02 

-.17266-02 

.12556-15 

-.80726-02 

.43226-04 

-.69046-15 
-, 12846-01 
.12496-02 

-.21106-14 

-  . 389 16  —  02 

.22916-02 


.44166-02 
•  «  62806-02 
. 79816-04 

-.38966-02 

. 49316-0? 
■  .  68 04t — 02 
-.32506-02 

.46816-02 

-.64646-02 

-.25096-02 

.90956-02 

-.13436-01 

-.30266-02 


.77436-02 

-.10426-01 

.43386-02 

•59396-02 

-.76356-02 

-.10546-02 

.99906-02 
-. 1 2666-01 
.12796-03 

.93526-03  - 

.04036-02 
-.71166-02 
.17146-02  - 


.62806-02 
.21 00E-03 
.  19*66-1* 

:${?? m 

.2008E-14 

.bBU*E-02 

.34526-02 

.13316-14 

.64646-02 

.26876-02 

-.62776-15 

-.12336-15 


10426-01 

44646-0? 

23856-14 

mm 

22606- 14 

18836-14 
288 7E-14 


.44746-02 

.62776-15 

-.79816-04 

.46706-02 

.16326-14 

•3B966-C2 

.47066-02 
-.11 J0t-14 
.32506-02 


.7o?0E-C? 

.1004£-»4 

-.43386-02 

-:fKff:ff 

-.1279E-03 

mm 


Iim:*S 

10596-14  .50266-02  -.17146-02 


Figure  30.  Partial  output  of  BOTZCW  for  Problem  2c. 
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Figure  31.  Parinl  output  of  BO'.TJVA  for  Prob  i  2c  when  the  loop  i)  added  to  the  capped  BOT. 
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Figure  32.  Partial  output  of  BOTRA  for  Problem  2c. 
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Figure  31’.  Concluded. 


2.6.3  Problem  3:  Wire  and  Aperture  Antennae  on  Wing  Section 

Consider  an  asymmetric  wing  section  shown  in  Figure  33a.  Two 
1/8X  monopoles  are  mounted  on  the  trailing  edge  of  the  wing  surface.  The 
monopole  at  the  center  of  the  surface  is  parasitic.  The  other  is  active.  The 
half-length  of  the  wing  is  1.38X. 

a)  Calculate  the  principal  plane  radiation  patterns  for  the  foregoing 
configuration.  Calculate  the  currents  on  the  BOT  surface. 

b)  Compute  the  electric  and  magnetic  near-fields  for  the  system  of 
radiators  at  selected  points  along  a  line  1/16X  from  the  wing  surface. 

c)  Calculate  the  coupling  between  the  two  monopoles  in  this  problem.  Next 
assume  that  there  are  three  aparture  antennas  (polarized  along  t)  on  the  wing 
surface,  centered  at  (tA,  zA)  and  (tfil  Zg)  and  having  an  axial  half-length  of 
La  and  Lg,  respectively.  Compute  the  coupling  between  these  slots. 

Solution  -  Figures  33b  and  33c  list  the  input  data  for  the  three  problems. 

The  open  BOT  system  matrix  is  calculated  by  executing  BOTZSS  and  BOTINV.  The 
wires  are  added  to  the  system  matrix  by  executing  BOTZ.iW,  BOTZVJW,  and 
B0T1NVA.  The  solution  to  Problems  3a  and  3b  is  obtained  by  executing  B0T«A 
(see  Figure  34a-34c  for  partial  outputs).  The  solution  to  Problem  3c  is 
obtained  by  executing  BOTAC  (see  Figures  35a-35b  for  partial  outputs). 
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Figure  34e.  Partial  output  of  BO  i  RA  for  Problem  3a. 
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Figure  36a.  Partial  output  of  BOTAC  for  Problem  3c. 
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Figure  35b,  Partial  output  of  BOTAC  for  Problem  3c. 


2.6.4  Problem  4:  Scattering  from  a  Closed  Cylinder 


Consider  a  circular  closed  cylinder  of  2.76X  length  and  0.216X  radius. 
(The  uncapped  cylinder  of  same  dimensions  Is  used  in  sample  Problem  1.) 

a)  Calculate  the  monostatic  (backscatter)  cross  section  of  the  body  for 
horizontal  (60)  and  veLtlcal  (<J>4>)  polarizations  as  a  function  of 
azimuthal  (6)  angles. 

b)  Calculate  the  bistatic  crois  section  of  the  body  for  horizontal  and 
vertical  polarization  as  a  function  of  0.  Assume  the  Illumination  is  normal 
(i.e.,  along  0  *  90). 

Solution  -  Figure  36  lists  the  input  data  for  the  two  problems.  The  closed 
BOT  system  matrix  is  calculated  by  executing  BOTZSC,  BOTZCC,  and  BOTINVA  using 
the  open  BOT  systems  matrix  from  Problem  1.  The  solution  to  Problem  4a  is 
obtained  by  executing  BOTSCM  (see  Figure  37  for  partial  output).  The  solution 
to  Problem  4b  is  obtained  by  executing  BOTSCB  (see  Figure  38  for  partial 
output) . 
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Figure  36.  Input  data  for  execution  of  Problems  4a4b. 
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Figure  37.  Partial  output  of  BOTSCM  for  Problem  4r.. 
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Figure  37.  Concluded. 
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Figure  38.  Partial  output  of  BOTSCB  for  Problem  4b. 
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3.  SYSTEMS  SECTION:  DETAILED  PROGRAM  DESCRIPTIONS 


A  detailed  description  of  each  A-STAR  program  (Figure  1)  is  given  below, 
including  a  description  of  the  flow  diagram,  subroutine  input/output 
arguments,  and  the  special  matrix  properties  used.  A  description  of  the 
common  variables  used  in  the  program,  along  with  storage  methods  used  in 
certain  arrays,  is  given  in  Appendix  A.  A  description  of  subroutines  is  given 
in  Section  3.6. 

3.1  BOTZSS  Program 

BOTZSS  generates  the  impedar.ee  submatrices  Z88^  for  modes  m,n  where 
0  <  m  <  NMODE  -  1,  -m  <  n  <  +  m,  and  |m  -  n|  <  NBAND.  The  impedance  matrices 
are  generated  in  the  lower  triangular  portion  of  ZgQT  irigure  39).  Symmetry 
conditions  are  then  applied  in  BOTINV  to  till  the  entire  ZgQT  matrix.  The 
structure  of  the  matrix  is  as  follows: 


9  8 

where  each  of  the  z  watrlces  la  comprised  of  four  subraatrlces  as  follows: 


-8Stt  1  „8StZ 

Z  |  Z 

m ,  a  i  ra,ri 

i 


SSZt  I  S9ZZ 

z  I  z 

m ,  n  |  m ,  n 


•  &Ww.V«nte  j. 


Figure  40  shows  the  flow  diagram  for  BOTZSS.  The  equation  numbers  refer 
to  the  analytical  expressions  in  Volume  I  of  this  report. 


The  computation  of  the  Green's  function  kernel  takes  advantage  of  the 
I 


G(i  +  ( j-1) j/2)  when  1  <  j 


G(j  +  (i-l)i/2)  when  i  >  j. 


Read  BOT  input,  compute 
parameters,  and  plot  the  BOT 
generating  curve 


Subroutines  BOTIN 
and  PLOTB 


Calculate  impedance  loading 
weights  for  the  BOT  and 
write  these  weights  to  disk 


Step  modes 
m  =  0  to  NMODR  -  1 
n  =  -  m  to  +  m 


Compute  Green's  kernel 
Gr.,n  and  Gmn  matri<:es 


Compute  impedance  matrix 
Z”n,  composed  of 

zsstt  ZSSU  zSSZt  d  Zsszz 
inn'  ‘mn  1  ‘mn  *  ana  £mn 


Write  Z^n  to  disk  and  line  printer 


X  All  \ 
modes 
complete 
\  / 


(Symmetry  conditions  applied) 


[Equations  (A-8),  (A-9),  (A-11),  <A-13)J* 


[Equations  (A  3)  to  ( A-S) )* 


'See  Volume  I 


Figure  40.  BOTZSS  flow  diagram. 
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The  section  that  computes  the  impedance  matrix  Zss  uses  the  following 

m,n 

symmetries : 


sstt 

m,n 


(zs8tt) 

tv  m,n 


» 


,  sstz . 

n  > 

t  m,n 


.  „sszt. 
-  n  (Z  ) 
m,n 


and 


2SSZZ  _  (z38Z2) 

m,n  t.  m,n 


Thus  only  the  upper  triangular  portion  of  each  of  the  Zs  ,  Z  ,  Z3  , 

8Szz  m,n  ®  ,r‘ 

and  Z  needs  to  be  computed.  The  remaining  portion  is  filled  us^ng  tie 
m,n 

symmetry  conditions  above. 

Since  the  z-directed  currents  on  the  BOT  are  expanded  as  t-2nTZ/L  _  n> 
the  n-0  mode  should  not  be  included.  However,  in  order  to  maiarain  a  parallel 
treatment  for  the  t-  and  z-directed  BOT  currents,  the  n-0  (z-directed;  mode  is 
included,  yielding 


sstz  sszt 

4.  =  Z 

mo  on 


sszz  sszz 

Z  -  Z  -  0. 
mo  ou 


This  addition  leads  to  a  singular  impedance  matrix  since  the  n-0  mode  z- 
dlrected  BOT  current  coefficients  are  arbitrary.  To  circumvent  this 
predicament,  the  BOTZSS  program  sets  Z33zz  -  I,  where  I  is  the  identity 
matrix,  and  the  matrix  inversion  can  be  performed  without  alteration,  which  is 
equivalent  to  forcing  the  n-0  (z-directed)  mode  BOT  current  coefficients  to  be 
zero. 


92 


3,2  BOTINV  Program 


* 


* 

* 


BOTINV  Ellis  the  ZB(yj,  matrix  using  the  output  file  from  B0T2SS  ar.d 
Inverts  this  matrix  according  to  the  user's  specifications.  The  following 
symmetries  are  <  sed  to  fill  from  the  partial  Zbq.j,  matrix  generated  by 

BOTZSS  (see  Figure  39  for  NMODE  =  4): 


Main  diagonal  symmetry 


(Compute  Z38  from  Zsa  ) 
n,m  m,n 


1  “ 

—  .  — 

zS8tt  j  z88t2 

zS8tt  (Z88Zt) 

n,ra  j  ntra 

a 

m,n  1  t  m,n 

1 

i 

_SSZt  1  „sszz 

1 

,_88tZ.|  _S8ZZ 

Z  1  Z 

-  (Z  )l  Z 

n,m  .  n,m 

t  ra,n  m,n 

—  1  _ 

-  1  __ 

Skeii  symmetry  (Compute  Z_n  _m  from  Z^^) 


7sstt 

-n,-m 

_ _ _ 

788tZ 

El 

sstt  .  .  sszt. 
m,n  ‘  t''  m,n  ’ 

I 

z88Zt 

788ZZ 

(Z8SC2)I  z88ZZ 
t  ro,n  .  m,n 

-n,-m 

-n,-a 

Flgi  re  41  shows  the  flow  diagram  for  BOTINV.  Three  types  of  matrix 
1*  ersiocs  are  allowed  In  BOTINV:  total  inversion,  main  diagonal  inversion, 
at 4  partla.  inversion.  Each  of  these  options  is  described  next. 

Total  Inve  slon  -  Total  inversion  is  performed  when  NBAND  >  2*NM0DE-1.  In 
this  case,  the  ZB0T  matrix  is  stored  by  columns  as  follows: 


(Z96tS  :  Z((n  +  NMODE- l)*LS2*(2*NMODE-l) 

m,n  ij 
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| 


t 


l 


* 


j 


-a**' 


Figure  4 1 .  BOTINV  flow  diagram 
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> 


+  ( j-l)*(2*NM0DE-l)*LS 


+  Cm  +  NMODE-1  )*LS  +  i) 


(Zsszt) 

is 

stored 

at 

<4"“) 

+  NM 

n>,n 

ra,n 

(zS8tz) 

is 

stored 

at 

(Z38tt) 

+  NM*LS*(2*NM0DE-1) 

«,n 

rafn 

ij 

,  sszzx 

.  sstzv 

(2  ) 

is 

stored 

at 

(Z  ) 

+  NM 

ra  „  n  ^  j 

nivn 

ij 

Once  the  ZfiQT 

matrix  is 

filled,  it 

is  inverted  using  Gaussian  elimination  with 

partial  pivoting,  and  written  to  disk  file  by  submatrices. 

Main  Diagonal  Inversion  -  Main  diagonal  inversion  is  performed  when  NBAND  -  1. 
In  this  case,  the  Individual  diagonal  submatrices  are  inverted  separately 
using  Gaussian  elimination.  For  m  *  0,  the  following  symmetry  is  used: 


tss 

— ra,-m 


— 

1  " 

i  ~ 

zsstt 

|  2SsCz 

zsstt 

■  _8BtZ 

1  ^  " 

-m ,  -m 

|  -ra,-m 

ra,m 

— 

•4*  — - 

m 

— 

1 - 

^sszt 

!  S8ZZ 

_88Zt 

|  _88ZZ 

L. 

2 

-z 

1  Z 

I  “lu  ,  *  iQ 

*a ,  in 

i  m ,  ro 

L 

1 

— 

i  _j 

which  implies  that 


ass  )‘ 

-m,-m 


f3S 

~m,-m 


has  the  same  symmetries.  Thus  only  z88  is  inverted  where  m  »  0  to  NMODE-1. 

F‘  ,m 

The  symmetries  are  used,  and  the  resulting  (2*NM0DE-1)  submatrices  are  written 
to  disk.  file. 
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Partial  Inversion  -  Partial  inversion  is  performed  when  NBAND<2*NM0DE-1 ,  and 


NBAND* 1.  In  this  case,  the  Z«aT  matrix  is  filled  only  with  the  Z 

oui  m,n 

submatrices  for  which  |m-nl  <  NBAND.  The  resulting  ZgQ^  matrix  has  a 

staircase-type  structure.  The  rest  of  the  matrix  is  sparse.  If  each  of  the 
s  s 

Zm  ^  submatrices  is  thought  of  as  an  individual  element,  Zg0«j,  can  be 
considered  as  a  banded  matrix.  A  modified  LU  (lower-upper  triangular) 
decomposition  can  then  be  used  with  all  arithmetic  operations  replaced  by  the 
corresponding  matrix  operations.  The  result  is  an  L  and  U  matrix  which  is 
also  of  a  staircase-type,  but  is  lower  and  upper  triangular,  respectively, 
when  the  submatrices  are  considered  as  individual  elements.  The  inverse  of 
ZBQT  can  be  computed  using  forward  and  backward  substitution,  again  replacing 
arithmetic  operations  with  matrix  operations.  The  result  is  a  full  inverted 
Zjjqt  matrix  which  is  written  to  disk  file  by  submatrices. 

8  S 

The  Zfl0T  matrix  is  stored  by  columns,  if  the  individual  submatrices  Z^  ^ 
are  considered  as  elements.  Only  the  banded  portion  is  stored.  When  NMODE  ■*  4 
and  NBAND  -  2  (refer  to  Figure  41),  Zg^^,  is  stored  in  the  following  order: 


Each  submatrix  is  stored  in  LS^  successive  locations  by  columns.  For  the 

ss  2 

example  above,  Z  ,  would  start  at  index  LS  +  1. 

9  J 


3.3  Description  of  Impedance  Generating  Programs 
3.3.1  BOTZSW  Program 

BOTZSW  generates  the  BOT -wire/ junction  impedance  submatrices  in  Figure  42 

for  -NMODE+1  <  m  <  NMODE-1.  Each  of  the  Zsw  matrices  is  comprised  of  up  to 

m 

four  submatrices  as  follows: 


z3w,t  I  zsj>t 

m  |  m 

_„T__ 

zsw*2  I  z8j’z 

m  ,  m 


where  t  and  z  correspond  to  the  t-  and  z-directed  cuirent  expansions  on  the 
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Figure  43  BOTZSW  flow  diagram. 


9S 


3.3.2  BOTZWW  Program 


BOTZWW  generates  the  wire/ junction-wire/ junction  impedance  submatrix  Z, 
which  is  comprised  of  up  to  four  submatricea  as  follows: 

zw  |  zwj 

» 

zjw  !  z^ 


where  w  ar.d  j  represent  the  wire  and  junction  expansions,  respectively. 

BOTZWW  usee  the  symmetry  relation  Z^w  »  tZw^,  w'lere  c  indicates  the  transpose 
operation.  Figure  44  shows  the  flow  diagram  for  BOTZWW. 


Read  wire/junction  input  SubrQutine  DATA|N 
and  compute  parameters 


s  Any 
junctions 
present 
v  i  ^ 


Compute  wire-junction  impedance 

matrix  and  us*  symmetry  to  obtain  (Table  1,  Volume  K 

junction  wire  impedance  elements 


Compute  junction-junction  n 

impedance  matrix 


[Table  1,  Volume  I] 


Compute  wire-wire  impedance 


Subroutine  ZLIST 


Figure  44.  BOTZWW  flow  diagram. 
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3.3.3  BOTZSC  Program 


SC 

BOTZSC  generates  the  BOT-cap/edge  impedance  submatrices  Z^  for 
-NMODE+1  <  m  <  NMODE-1.  Each  of  the  Zm  matrices  is  comprised  of  up  to  six 
submatrices  as  follows: 


zsc,tt  I  2sc,tp  .  2SC»te 

m  |  m  I  m 

- h - 4- - . 

Zsc,zt  |  Zsc,zp  j  Zsc,ze 
m  |  m  |  m 

where  the  third  superscript  t  or  z  T°f_ca  to  the  t-  or  z-directed  expansions 
on  the  BOT,  respectively;  the  fourth  superscript  refers  to  the  t  or  p 
expansions  on  the  cap,  and  e  refers  to  the  edge  expansions.  Figure  45  shows 
the  flow  diagram  for  BOTZSC. 

3.3.4  BOTZCC  Program 

BOTZCC  generates  the  cap/edge-cap/edge  impedance  submatrix  Zcc,  which  is 
comprised  of  up  to  nine  submaurlces  as  follows: 


— * 

“ 

zCC,tt 

zcc,tp 

zcc,te 

“cTIpt 

2cc,pp 

zcc,pe 

zcc,et 

cc,ep 

2 

_cc.ee 

r 


where  the  t,  p,  and  e  represent  t,  p  cap,  and  edge  expansions,  respectively. 
BOTZCC  uses  the  symmetry  relations  Zcc,et  -  cZcc*te  and  Zcc,ep  -  t zcc,pe , 
where  t  indicates  the  transpose  operation.  Figure  46  shows  the  flow  diagram 
for  BOTZCC. 
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Read  Btf  input,  compute 
parameteis,  and  plot  the  BOT 
geneiating  curve 


Subroutines  BOTIN 
and  PLOTS 


Read  cap/edge  input  end 
compute  parameters 


Subroutine  CAPIN 


Step  modes  m  =  —  NMGDE  +  1 
to  NMODE  -  1 


Compute  Green's  kernel  matrix 
•  between  BOT  and  cap/edgo  patches 


Compute  BOT-cap  impedance 


[Equation  (A-20a)  -  (A-21)l* 


^  Cap 
edges 
present 
\  ?  / 


Compute  BOT-cap  impedance 
matrix 


[Equations  (A-31a)  —  (A-31b))“ 


-  All  > 
modes 
complete 
\  / 


“See  Volume  I 


Write  to  disk  and  line  printer  Subroutine  ZLIST 


Figure  45.  BOTZSC  flow  diagram. 
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[ 


Read  BOT  input,  compute  Subroutines  BOTIN 
parameters,  and  plot  BOT  and  plOTB 
generating  curve  _ 


I 


Read  cap/edge  input  and  |  Subroutine  CAPIN 
compute  parameters 


Calculate  surface  impedance 
loading  weights  for  cap/edge 
and  write  these  weights  to 
disk 


ComputeTireen's  kernel  self-term  [Equati0n  (A-17>1‘ 
matrix  GS  for  cap  P«cl«s _ 

rx- ^ 

Z  I  *  I 


I  matrices 

Wccpp 


/  Are 
/  edge 
No  /  transition 

terms 

included 

7 


-edge  impedanc  .  1  [Equations  (A-32a)  (A  32b)] 

and  Z** ^ _ J 

y  matrix  Zee - — — * 


f~  Compute  cap 
1  submatnces  Z  1  and  L 

I - 


U|  Write  Zcc  to  disk  and  line  Subroutine 

|  printer  _ _  - 


ZLIST 


•See  Volume  I 


STOP 


Figure  46.  BQTZCC  flow  diagram. 
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3.3.5  BOTZCW  Program 


BOTZCW  generates 
is  comprised  of  up  to 


the  cap/edge -wire/ Junction  impedance  submatrLx  Z,  which 
six  subm&trices  as  follows: 


zcw,p  |  zcj,p 
ew  T~  e  j 

‘  I  * 


where  c,  p,  and  e  refer  to  the  t,  p  cap,  and  edge  expansions,  respectively, 
and  w  and  j  refer  to  the  wire  and  junction  expansions,  respectively.  The 
interaction  between  the  cap/edge  currents  and  the  junction  currents  is  set  to 
zero  (i.e.,  zc^,t:  =  Zc^,p  ■*  =>  0).  Figure  47  shows  the  flow  diagram  for 

BOTZCW. 


3.4  BOT1NVA  Program 

B0T1NVA  Inverts  the  system  matrix  when  either  wires  or  caps  are  added  to 
a  BOT  using  a  previously  inverted  system  matrix.  The  new  system  matrix  to  be 
inverted  can  be  written  in  partitioned  form  as: 


where  P  is  the  old  system  matrix  for  which  P~*  already  exists  on  disk,  and 
where  Q,  R,  and  S  are  impedance  matrices  arising  from  the  addition  of  wires 
and/or  caps.  The  inverse  of  this  system  is  of  the  form: 


i 


Figur*  47.  BOTZCW  flow  diaptm. 


e 


l>JO  X  junctions 
"  \  present 

X  ?  ^ 


Set  Zej  =  0 


Write  Z  to  disk  and  line 
printer 


Subroutine  ZUST 


where 


N  -  (S-R?'*lQ)-1 
L  -  ~P_1QN 
M  -  -NRP’1 
K  -  P_1(I-QM). 


Figure  47.  Concluded. 


Since  P~1  already  exists,  the  new  inverse  can  be  found  by  matrix  multipli¬ 
cation  and  by  calculating  the  inverse  of  a  matrix  having  the  same  order  as 
S.  Figure  39  shows  the  system  matrix  that  is  obtained  when  caps  are  added  to 
an  old  system  matrix  containing  an  opan  BOT  with  wires  and  Junctions.  The  S 
matrix  corresponds  to  the  impedance  matrix  generated  by  BOTZCC  for  this 
case.  The  same  system  could  be  generated  by  adding  the  caps  first,  which 
woula  result  in  a  rearrangement  of  the  matrices.  In  either  case,  BOTINVA 
arranges  the  matrices  properly.  Figure  48  shows  the  flow  diagram  for  BOTINVA 


I. 


<•> 

0 


Figure  48.  BOTINVA  flow  diagram. 
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Subroutine  1NVPAR 


Subroutine  PUTPQRS 


Fiii«r«s8,  Concluded. 


3.5  Radiation/ Scattering  Analysis  Programs 

BOTRA  Program  -  BOTRA  computes  the  radiated  far  and  near  fields  resulting 
from  slot  and/or  wire  antennas  on  the  BOT  with  or  without  end-caps  present 


(Sections  5  and  6,  Volume  I).  The  presence  of  wires  or  caps  in  the  inverted 
system  matrix  file  is  indicated  by  assigning  pseudo  mode  numbers  to  the  wire 
and  cap  blocks  in  the  matrix  (see  Figure  39).  Figure  49  shows  the  flow 
diagram  for  BOTRA. 


BOTSGB  Program  -  BOTSCB  computes  scattered  far  and  near  fields  in  the 
bistatic  mode,  with  or  without  wires  and  arid-caps  present  (Section  5.3,  Volume 
I).  Figure  50  shows  tne  flow  diagram  for  BOTSCB. 


BOTSCM  Program  -  BOTSCM  computes  scattered  far  fields  in  the  monostatic 
mode  (Section  5.3,  Volume  I).  Figure  51  shows  the  flow  diagram  for  BOTSCM. 
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BOTAC  Program  -BOTAC  computes  the  antenna  coupling  between 
antennas  located  on  th<  BOT  surface  (Section  6.3,  Volume  I). 
the  flow  diagram  for  BOTAC. 


Figure  49.  BOTRA  flow  diagram. 


slot  and  wire 
Figure  52  ahr 
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Subroutine  MUL.TVC 


Subroutines  RBOT, 
RWIRE.  and  RCAP 
(Table  2.  Volume  I] 


Subroutine  MULTYR 
[Equation  (A- 34), 
Volume  I) 


Figure  49.  Continued. 


Figure  49.  Concluded. 
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Figure  50.  BOTSCB  flow  diagram. 


t 

i 


A 


£ 


DO  90  IANG  =  1.  NANG 

1 

Set  the  and  6S  scattering  angles 

DO  110K«  1.  NT ;  fix  4>i  and  1 

Calculate  transfer  matrices 
for  $  and 

T _ 

Ca'culate  2  Rm(ds.  0$)  I,^  for  66. 
ip4>.  8<t>,  and  $0  polarizations 

1 

Print  RCS  for  <t>%  and 

; 

110  CONTINUE 

1 

90  CONTINUE 

1 

List  and  plot  currents 

I 


Oo  490  I  TEST  *  1.  NTEST 

I 


Compute  modal  current 
coefficients  (for  x,  y,  and  z 
components)  for  the  electic 
and  magnetic  near-fields. 
Sum  ewer  all  modes  to 
calculate  the  electric-  and 
magnetic-field  components 
resulting  from  6  and  <p 
polarized  incident  waves 

j  Print  near-field  results 

I 

490  CONTINUE 


STOP 


Subroutines  RBOT, 
R WIRE,  and  RCAP 
[Table  2,  Volume II 


Subroutine 

MULTYR 


Subroutines  LPCUR 
and  PLOTC 


(Near-field  analysis) 


Subroutines  NEARB, 
NEARC,  and 
NEARW 

[Equations  47,48.50, 
52,55,  and  o6'j* 


•See  Volume  I 


Figure  50.  Concluded. 
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(Equation  4  J]* 

Subroutines  MULTYR, 
and  MULTRC 


•See  Volume  I 


Figure  SI.  BOTSCM  flow  diagram. 
113 


Figure  62.  BOTAC  flow  diagram. 


3.6  Subroutine  Descriptions 

A  description  of  the  subroutines,  called  by  each  of  the  A-STAR  ptograms 
(Figure  1),  follows.  These  subroutines  are  referenced  in  the  flow  diagrams 
for  the  program  descriptions.  The  variable  descriptions  used  In  these 
routines  are  given  in  Appendix  A. 

BOTIN  -  Subroutine  BOTIN  reads  the  BOT  input  geometry  from  the  user  input  data 
(see  Section  2.3.1),  plots  the  BOT  generating  curve,  and  computes  all  BOT 
segmentation  parameters  used  in  calculations  involving  the  BOT. 

CAP IN  -  Subrout in;  CAPIN  reads  the  cap/edge  input  geometry  from  the  user  input 
data  (see  Section  2.3.2)  and  computes  all  cap  and  edge  parameters  used  in 
calculations  involving  the  caps  and  edges. 


CSIMP  -  Subroutine  CSLMP  is  a  Simpson  Integration  routine  with  a  calling 


statement : 


CALL  CSIMP(F  ,A,  8  ,DEL,  IMAX,  SI  1 ,  S ,N, IER)  . 

T» 

This  routine  computes  S  =>  F(x)dx  using  the  method  of  successive  bisections 

of  the  interval  until  either  a  relative  error  or  DEL  is  achieved  or  IMAX 
bisections  have  been  performed.  F  must  be  declared  external  in  the  calling 
program.  The  following  are  returned  by  CSIMP: 


S  -  Approximate  value  of  the  integral. 

Sll  -  Previous  approximation  to  the  Integral.  Convergence  has  occurred  if 


S  -  SI1 
S 


<  DEL. 


N  -  Number  of  Intervals  used  in  computing  S. 

IER  -  Error  return.  IER  -  0  indicates  that  convergence  has  occurred. 


DATAIN  -  Subroutine  DATA1N  skips  over  the  BOT  and  cap/edge  input  data,  reads 
the  wire/ junction  input  geometry  from  the  user  input  data  (see  Section  2.3.3), 
and  computes  all  wire  and  junction  segmentation  parameters  used  in 
calculations  involving  the  wires  and  junctions. 

GETPQRS  -  Subroutine  GETPQRS  retrieves  the  partitioned  system  matrix^  ^jfrom 
disk,  using  the  following  unit  numbers: 

P  1  is  read  on  unit  1 
S  is  read  on  unit  4 

Q  13  read  on  unit  2,  and  if  both  wires  and  caps  are  present  In  the  system 
matrix,  an  additional  part  of  Q  is  read  on  unit  3. 


The  matrix  R  is  obtained  as  the  transpose  of  Q. 


GREEN  -  Subroutine  GREEN  calculates  the  Green's  function  kernel  used  in  the 
calculation  of  wire-wire  impedance  matrix  elements. 


■tt 


4 

i 

i 

t 

l 
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INVBAN  -  Subroutine  INVBAN  is  a  modification  of  a  standard  banded  matrix 


inversion  routine  using  LU  decomposition  without  pivoting,  where  only  the 
banded  portion  is  stored  by  columns.  Arithmetic  operations  were  replaced  by 
their  corresponding  matrix  operations,  and  indices  were  multiplied  by  LS2 
since  the  elements  were  replaced  by  matrices.  The  calling  statement  and 
arguments  follow: 


CALL  INVBAN(LS , NMODE , NBAND , NZ , A, Z , WORK) , 

where  LS,  NMODE  and  NBAND  are  described  in  Appendix  A. 

INPUT:  NZ  -  Array  used  for  indexing.  In  a  normal  banded  matrix  A, 

NZ(I)  *■  NZ(l-l)  +  (the  number  of  zeros  below  the  band  in 
column  1-1)  +  (the  number  of  zeros  above  the  band  in  column 
I),  where  NZ(1)  «  0.  If  A  has  order  i  with  only  the  banded 
portion  stored  by  columns,  then  A^  wil^  be  stored  in  location 
n(j-l)  +  i  -  NZ(j). 

A  -  Array  containing  the  staircase-type  matrix  to  be  inverted, 
with  storage  details  described  above. 

Z  -  Array  of  length  LS2  used  as  work  area. 

WORK  -  Array  of  length  LS  used  as  work  area. 

In  addition,  three  variables  are  passed  in  common  as  follows: 

COMMON  NM,JK(4),LR, 
where 

NM  -  Number  of  triangle  functions 

JK  -  Work  array  of  length  4 

LR  -  Work  array  of  length  LS. 

INVPAR  -  Subroutine  lN^PAR  inverts  a  partitioned  matrix  of  the  form  i  a- 
where  P  *  already  exists.  The  calling  statement  and  arguments  follow: 

CALL  1NVPAR(P1,Q,R,S,W,LR,N,M). 

INPUT:  PI  -  Matrix  containing  P~*,  and  on  return  it  contains  the 
partitioned  part  of  the  Inverse. 

Q,R,S,  -  Submatrices  in  the  partitioned  matrix,  and  on  return  they 
contain  the  submatrices  of  the  Inverse. 
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W  -  Complex  work  array  of  length  max  (N,M). 

LR  -  Work  array  of  length  H. 

N  -  Order  or  matrix  P. 

M  -  Order  of  matrix  S. 

LINEQ  -  Subroutine  LINEQ  Is  a  standard  matrix  Inversion  routine  using  Gaussian 
elimination  with  partial  pivoting,  with  the  following  calling  statement  and 
arguments: 

CALL  LINEQ (LL.C.LR), 

where 

LL  -  Order  of  matrix  to  be  Inverted. 

C  -  Array  containing  the  matrix  to  be  inverted,  stored  by  columns.  On 
output,  C  contains  the  inverted  matrix. 

LR  -  Array  of  length  LL  used  as  a  work  space  during  the  pivoting  process. 

LIST  -  Subroutine  LIST  prints  Individual  Yffl  n  submatrices  cn  the  line  printer 
and  writes  them  to  a  disk  file. 

LISTA  -  Subroutine  LISTA  prints  individual  Y^  n  submatrices  that  correspond  to 
wires  and/or  caps,  and  writes  them  to  a  disk  file. 

LPCUR  -  Subroutine  LPCUR  lists  and  plots  the  BOT  currents  (magnitude  and  phase 
as  a  function  of  Z/BL)  for  each  triangle  function  on  the  body,  and  then  lists 
cap/edge  and  wire/ junction  currents. 

NEARB  -  Subroutine  NEARB  calculates  the  electric  and  magnetic  current 
coefficients  for  one  mode  of  the  BOT  current  expansion,  and  for  one  near-field 
test  point. 

NEARC  -  Subroutine  NEARC  calculates  the  electric  and  magnetic  current 
coefficients  for  the  caps  resulting  from  one  near-fleld  test  point. 

NEARW  -  Subroutine  NEARW  calculates  the  electric  and  magnetic  current 
coefficients  for  the  wires  resulting  from  one  near-field  test  point. 
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ORDER  -  Subroutine  ORDER  checks  that  the  user  input  data  agree  with  the 
parameters  contained  in  the  system  matrix  disk  file,  and  prints  the  order  in 
which  caps  and/or  wires  were  added  to  the  system. 

PLOT  -  Subroutine  PLOT  generates  an  x  vs.  y  plot  on  the  line  printer. 

PLOTB  —  Subroutine  PLOTB  plots  the  points  on  the  generating  curve  of  the 
/iOT.  Points  on  the  BOT  corresponding  to  triangle  function  peaks  are  indicated 
with  a  plus  sign.  The  calling  statement  is 

CALL  PLOTB(X,Y,N,NR), 

where 

X  -  Array  of  x  coordinates  to  be  plotted 
Y  -  Array  of  y  coordinates  to  be  plotted 
N  -  Number  of  points  to  be  plotted 

NR  -  Number  of  line  printer  rows  to  use  for  the  y-axis. 

The  routine  uses  51  columns  for  the  x-axis,  with  the  dynamic  range  on  both  the 
x  and  y  axes  equal.  Hence,  depending  upon  the  type  of  line  printer  used,  NR 
may  have  to  be  adjusted  to  obtain  a  plot  that  is  not  distorted  (i.e.»  the  x 
and  y  axes  have  approximately  the  same  physical  length  on  the  line  printer 
output). 


PLOTC  -  Subroutine  PLOTC  plots  the  magnitude  and  phase  of  the  currents  on  a 
given  BOT  triangle  function.  The  calling  statement  is  as  follows i 


where 


CALL  FL0TC(Y1,Y2), 


IXrOT:  Y;  -  Array  containing  the  current  magnitude  at  41  equally  spaced 
z  coordinates. 

\l  ~  Array  containing  the  current  phase  at  41  equally  spaced 
z  coordinates. 


PUT&I-S  “  Subroutine  PUTPQRS  writes  the  partitioned  inverted  system  matrix  to 

a  dlt/k  file  by  submatrices  Ym  „  where  m  and  n  are  either  BOT  mode  numbers  or 

m,u, 

pseudo  mcda  numbers  corresponding  to  wires  and  rep«. 
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RBOT  -  Subroutine  RBOT  computes  the  mode-independent  part  of  the  BOT  transfer 
matrices  for  ($,8)  angles. 

RCAP  -  Subroutine  RCAP  computes  the  cap/edge  transfer  matrices  for  C$,6) 
angles. 

RWIRE  -  Subroutine  RWIRE  computes  the  wire/ Junction  transfer  matrices  for 
(4>,8)  angles. 

SB0T1N  -  Subroutine  SBOTIN  skips  over  the  BOT  geometry  input  data. 

S CAP IN  -  Subroutine  SCAPIN  skips  over  the  cap/edge  input  data. 

SLOTIN  -  Subroutine  SbOTIN  reads  the  slot  antenna  input  data  (see  Section 

2.3.4). 

VBOT  -  Subroutine  VBOT  computes  the  BOT  voltage  array  for  a  given  mode  number 
resulting  from  the  slot  antennas  on  the  BOT. 

VWIRK  -  Subroutine  VWIRE  reads  the  wire/ junction  voltages  and  generates  the 
wire/ junction  voltage  array. 

VIR  IN  -  Subroutine  WIREIN  reads  the  wire/ junction  input  geometry  from  the 
user  Input  dato  (see  Section  2.3.3)  and  computes  all  wire  and  junction 
segmentation  parameters  used  in  calculations  involving  the  wires  and 
junctions. 

ZL1ST  -  Subroutine  ZLIST  prints  impedance  matrices  on  the  line  printer  by 
submatrices  according  to  the  type  of  current  expansions  contained  in  the 
impedance  matrix. 
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APPENDIX  A:  DICTIONARY  OF  COMMON  PROGRAM  VARIABLES 
(Input  variables  are  described  In  Section  2.3;  equation  numbers  refer 
to  expressions  in  Volume  I;  page  numbers  refer  to  Volume  II) 


AC(I) 


ANC 

ANG1 

ANC2 

BK 

BKL 

BL 

CB 


CBP 


CBT 


CC 

CCP 

CCT 

CPC(I) 

CV(I) 

cw 


-  Area  of  cap  triangle  formed  by  connecting  the  BOT  segment  I  with 
the  cap  center  point  (XC,  YC).  The  area  of  the  Jth  trapezoidal 
patch  on  this  cap  triangle  is  then  given  by  AC(I)*(RH0C(J+1)**2 
-RHOC(J)**2) . 

-  Input  array  of  fixed  radiation  or  scattering  angles  (page  15). 

-  Input  array  of  starting  radiation  or  scattering  angles  (page  15). 

-  Input  array  of  ending  radiation  or  scattering  angles  (page  16). 

-  Wavenumber  (meters-*). 

-  BK*BL. 

-  Half  length  of  BOT  (meters). 

-  Array  containing  modal  t-  and  z-d5.rected  currents  on  the  BOT 
[Equation  (2)]«  CB  [(m  +  NMODE  -  1)*LS  +  J]  contains  the  t- 
directed  current  for  mode  m  on  triangle  function  J.  The  z-dlrected 
current  is  stored  in  CB  [(m  +  NMODE  -1)*LS  +  J+NM].  Used  in  BOTRA. 

-  Array  containing  modal  t-  and  z-directed  currents  on  the  BOT  for  a 
^-polarized  incident  wave.  See  variable  CB  for  storage  details. 
Used  in  BOTSCB  and  BOTSCM. 


-  Array  containing  modal  t-  and  z-directed  currents  on  the  BOT  for  a 
6-polarized  incident  wave.  See  variable  CB  for  storage  details. 
Used  in  BOTSCB  and  BOTSCM. 


-  Array  containing  t-  and  p-directed  currents  on  the  caps 
currents. 

-  Array  containing  t-  and  p-directed  currents  on  the  caps 
currents  for  a  ^-polarized  incident  wave. 

-  Array  containing  t-  and  p-directed  currents  on  the  caps 
currents  for  a  6-polarized  incident  wave. 

-  Cosine  of  angle  for  BOT  generating  curve  segment  1. 

3  of  Volume  I.)  Used  when  caps  are  present. 


2  of  Volume  1.)  CV(J)  corresponds  to  v  on  segment  J. 
-  Array  containing  wire  and  Junction  currents. 


and  edge 
and  edge 
and  edge 
(See  Figure 
(See  Figure 
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CWP 

CWT 

DH(I) 

DHW(I) 

DRHOC(I) 

DTOR 

DXW(I) 

DYW(I) 

DZW(I) 

ESC 


ESCP 


ESCT 


EO 

ETA 

EWGHT 

G 


GB 

GP 

GS 

GT 


Array  containing  wire  and  junction  currents  for  a  ^-polarized 
In' i dent  wave. 

Array  containing  wire  and  junction  currents  for  a  9-polarlzed 
Incident  wave. 

Length  of  generating  curve  segment  I  (meters). 

Length  of  wire  segment  I  (meters). 

Length  of  normalized  radial  segment  I  on  the  caps. 
n/l80°. 

x  coordinate  variation  for  wire  segment  1  (meters), 

y  coordinate  variation  for  wire  segment  I  (meters), 

z  coordinate  variation  for  wire  segment  I  (meters). 

Array  containing  electric  near-field  radiation  components  of 
£(r')  in  the  x,  y,  and  z  directions  [Equations  (46-50)]  stored  In 
ESC(l-3),  respectively.  Used  In  30TRA. 

Array  containing  electric  near-field  scattering  components  result¬ 
ing  from  a  ^-polarized  incident  wave  [Equations  (46-50)].  Used  In 
BOTSCB . 

-  Array  containing  electric  near-field  scattering  components  result¬ 
ing  from  a  0-polarlzed  incident  wave  [Equations  (46-50)].  Used  In 
BOTSCB. 

-  Input  array  for  slot  antenna  excitation  [Equation  (39)]  (page  14). 

-  n  ■  7e  «  376.707  i). 

o  o 

-  Array  of  weights  used  for  impedance  loads  on  the  edges. 

-  Array  containing  the  Integrated  Green's  function  kernel  G,^ 
[Equations  (A-8)  to  A-12)].  In  BOTZSS,  G  is  symmetric  with  only 
the  upper  triangular  portion  stored  by  columns  from  Index  1  to 
(NP  -  l)*NP/2. 

G^j  is  stored  In  location  G[i  +  (J  -  l)*j/2]  when  i  <  j. 

-  Array  containing  the  Integrated  Green's  function  kernel  G^ 
[Equation  (A-6a)].  See  variable  G  for  storage  details. 

-  Array  containing  the  fields  for  the  NT  radiation  angles 
to  compute  4>-polarlzed  gain. 

-  Array  containing  the  Integrated  Green's  function  kernel  self -terms 
on  the  caps. 

-  Array  containing  the  fields  for  the  NT  radiation  angles  to 
compute  0-polarlzed  gain. 
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Gl  -  Array  containing  the  integrated  Green's  function  kernel  G 

m 

for  caps  and  wires. 

G2  -  Array  containing  the  integrated  Green' b  function  kernel  G 

m 

for  caps  and  wires. 

GIE  -  Arvay  containing  the  Integrated  Green's  function  kernel  G 

in 

for  edges. 

G2E  -  Array  containing  the  integrated  Green's  function  kernel  G 

m 

for  edges. 

HO  -  Array  containing  the  Integrated  Green's  function  kernel 

for  the  magnetic,  near  fields  [Equation  (52)]. 

HI  -  Array  containing  the  integrated  Green's  function  kernel 

for  the  magnetic  near  fields  [Equation  (52) j. 

HSC  -  Array  containing  magnetic  near-field  radiation  components  in  the 

x,  y,  and  z  directions  [Equations  (52)— (56) ] .  HSC  •  S(r'),  used  in 
BOTRA. 

HSCP  -  Array  containing  magnetic  near-field  scattering  components 
resulting  from  a  ^-polarized  incident  wave.  HSCP  *  H(r'), 
used  in  BOTSCB. 

HSCT  -  Array  containing  magnetic  near-field  scattering  components  result¬ 
ing  from  a  9-polarlzed  incident  wave.  HSCT  -  H(r'), 
used  in  BOTSCB. 

I EDGE  -  Indicates  whether  the  BCT  generating  curve  is  open  or  closed. 

IEDGE  *  0  for  closed  and  1  for  open. 

INDJW  -  Input  array  indicating  which  points  in  the  wire  array  are  junction 
points  (page  12). 

1NDTJ(I)  -  The  wire  segment  index  at  which  the  I-th  junction  half  triangle 
function  starts. 

INDTW(I)  -  The  wire  segment  index  at  which  the  I-th  wir^  triangle  function 
starts. 

INDW  -  Input  array  indicating  the  wire  indices  at  which  each  of  the  NW 
wires  start  (page  12). 

1PCANE  -  Input  array  Indicating  whether  corresponding  element  of  array  ANG 
Is  a  0  or  $  angle  (page  15). 

IS  -  Input  array  for  specifying  location  of  slot  antennas.  [See 

Equations  (36)-(39).]  (page  13). 
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IUNIF  -  Indicate  whether  the  BOT  generating  curve  has  uniform  or  nonuniform 
segmentation, 

IUNIF  -  1  for  uniform  segmentation  and  IUNIF  »  0  tor  nonuniform 
segmentation. 

KG  -  NP  -  1. 

LC  -  Total  number  of  two-dimensional  triangle  peaks  on  the  caps. 

LC2  -  2*LC. 

LE  -  Total  number  of  triangle  peaks  on  the  edges. 

LR  -  (NPR-3)/2  *  the  number  of  triangle  functions  on  one  of  the  caps  in 

the  radial  direction. 

LS  -  Order  of  each  submatrix.  LS  -  NP  -  3. 

LSS  -  LS*LS. 

I.W  -  (NPW-3*NW)/2  ■  the  total  number  of  triangle  functions  on  the  wires. 

LWJ  -  LW+NJ. 

M  -  Mode  number  m. 

MORD  -  Two-dimensional  array  that  indicates  the  order  in  which  wires  and 
ca,s  have  been  added  to  the  system  matrix.  MORD(1)-0  if  caps  are 
not  contained  in  the  system  matrix.  If  caps  are  contained  in  the 
system,  MORD(l)  .ontains  the  pseudo  mode  number  corresponding  to 
the  location  of  the  caps  in  the  system  matrix  (i.e.,  MORD( 1)-NM0DE 
if  caps  were  added  first,  and  MORD( 1)-NM0DE+1  if  caps  were  added 
second).  Similarly,  M0RD(2)-0  if  wires  are  not  contained  in  the 
system  matrix,  and  if  wires  are  contained  in  the  system  matrix, 
M0RD(2)  contains  the  psuedo  node  number  for  the  wires. 

N  -  Mode  number  n. 

NANG  -  Input  variable.  Number  of  fixed  radiation  or  scattering  angles 
(page  15). 

NBAND  -  Input  variable.  Number  of  submatrix  diagonal  bands  used  in  Z  JL 

BUT 

(page  9). 

NC  -  Number  of  caps. 

NE  -  Indicates  whether  cap/edge  terms  are  included.  NE“0  if  edge  terms 

are  not  present,  and  NE-NC  if  edge  terms  are  present. 

NJ  -  Number  of  junctions. 

NM  -  Number  of  triangle  functions  on  the  BOT  generating  curve. 

NM  -  (NP-3)/2, 

NMODE  -  Input  variable.  Number  of  nonnegative  modes  (page  9). 
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1 


NM2 

NM4 

NP 

NPR 

NSP 

NPW 

NSA 

NX 

NTESX 

NW 

PHII 

PHIR(K) 

PHIS(K) 

RADD(I) 

RADJ(I) 

RADW(I) 

RBP 


RBT 


-  Qtder  of  each  submatrix.  NM2-NP-3. 

-  NM*4. 

-  Number  of  points  on  the  BOX  generating  curve.  (See  Section  2.3.1.) 

-  Number  of  radial  points  on  each  cap. 

-  Input  variable.  Number  of  diagonal  bands  used  in  BOX  impedance 
matrices  (page  9). 

-  Xotal  number  of  points  on  the  wires. 

-  Number  of  slot  antennas  on  the  BOT. 

-  Input  variable.  Number  of  radiation  and  scattering  angles  (page 
15). 

-  Input  variable.  Number  of  test  points  for  near-flelds  (page  16). 

-  Number  of  wires. 

-  4>  angle  for  the  Incident  wave  (degrees).  PHII  *  In  Equation 
(43). 

-  <f>  angle  for  the  radiated  fields  (degrees).  PHIR(K)  »  <t>  in 
Equation  (32). 

-  <j>  angle  for  the  scattered  fields  (degrees).  PHIS(K)  =*  in 

s 

Equation  (43). 

-  Radius  of  the  I-th  junction  disk  (meters). 

~  Radius  of  the  I-th  junction  wire  (meters). 

-  Radius  of  I-th  wire  (meters). 

-  Array  containing  the  mode-independent  portion  of  the  BOX  R^  and 
Rz^  matrices  (i.e.,  with  the  a  term  removed)  resulting  from  a  <j>- 
polarized  incident  wave.  The  order  of  storage  Is 

KV * -  1  to  NM}  followed  by 
{(R^)l.  1  -  1  to  NM}  . 

RBP  can  contain  the  transfer  matrices  for  several  ($,  6)  angles. 

In  this  case,  the  starting  index  is  offset  by  a  multiple  of  2*NM. 

t  0 

-  Array  containing  the  mode-independent  portion  of  the  BOX  R^  and 

RZ®  matrices  (i.e.,  with  the  a  term  removed),  resulting  from  a 
n 

6-polarlzed  incident  wave.  Xhe  order  of  storage  Is 

{(R^0)^,  1  =  1  to  NM}  followed  by 

{ (Rz0)  ,  i  «  1  to  NM}  . 
n  i 

As  in  RBP,  the  starting  index  can  be  offset  by  a  multiple  of 
2*NM. 
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RC(l) 

RC1(I) 

RCP 


RCT 

RHOC 

RHOC1 

RWP 


RWT 

RWGHT 

SPC(I) 

SPP 

SPT 

STP 

STT 

SV(I) 


T 


-  RadtaL  distance  from  the  center  of  the  cap  to  the  1-th  BOT  point 
(meters) . 

-  Radial  distance  from  the  center  of  the  cap  to  the  I-th  BOT  segment 
(meters) . 

-  Array  containing  the  cap/edge  transfer  matrices  resulting  from  a 
^-polarized  incident  wave.  The  order  of  storage  is 

{R^*,  i  *  1  to  LC}  followed  by 
{R^,  i  -  1  to  IX }  followed  by 
{R®*,  1  =*  1,  LE}. 

RCP  can  contain  the  transfer  matrices  for  several  (<f>,  6)  angles. 

In  this  case,  the  starring  index  is  offset  by  a  multiple  of  2*LC+LE. 

-  Array  containing  the  cap/edge  transfer  matrices  resulting  from  a 
0-polarized  incident  wave.  See  variable  RCP  for  storage  details. 

-  Input  array  of  normalized  radial  coordinates  on  each  cap. 

-  Normalized  radial  coordinates  for  the  radial  segments  on  each  cap. 

-  Array  containing  the  wire/ junction  transfer  matrices  resulting  from 
a  ijp-polarized  incident  wave.  The  order  of  storage  is 

{RW*,  i  **  1  to  LW}  followed  by 
{R^*,  i  -  1  to  NJ(. 

RWP  can  contain  the  transfer  matrices  for  several  (<)>,  0)  angles. 

In  this  case,  the  starting  index  is  offset  by  a  multiple  of  LW+NJ. 

-  Array  containing  the  wire/ junction  transfer  matrices  resulting  from 
a  9-polarized  incident  wave.  See  variable  RWP  for  storage  details. 

-  Array  of  weights  to  be  used  for  p-directed  Impedance  loads  on  the 
caps. 

-  Sine  of  *  angle  for  BOT  generating  curve  segment  I.  (See  Figure  3 

P 

of  Volume  I).  Used  when  caps  are  present. 

-  Array  containing  o1^  for  the  NT  scattering  angles  [Equation  43)]. 

-  Array  containing  for  the  NT  scattering  angles  [Equation  43)]. 

-  Array  containing  for  the  NT  scattering  angles  [Equation  (43)]. 

00 

-  Array  containing  o  for  the  NT  scattering  angles  [Equation  (43) J. 


-  Sine  of  angle  for  BOT  generating  curve  segment  I.  (See  Figure  2 

of  Volume  I.)  [Equation  (A-3).]  SV(J)  corresponds  to  v  on 

P 

stgraent.  J. 


-  Array  containing  the  values  of  the  triangle  functions  TC. 

P 

T[(K  -  1)*4  +  p]  contains  the  value  of  the  k-th  triangle  function 


TBE 

TCE 

TCR 

TCT 

TEXC 
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THS(K) 

TJ 

TP 

TPBE 

TPCE 
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TPCT 

TPJ 

TPW 

TW 

TVJGHT 


over  the  p-th  segment  forming  It.  1  <  p  <  4.  [Figure  2,  Volume 

I]- 

-  Array  con.aining  the  values  of  the  edge  half-triangle  functions  on 
the  BOT. 

-  Array  containing  the  values  of  the  edge  half-triangle  functions  on 
the  cap. 

-  Array  containing  the  values  of  the  radially  directed  triangle 
functions  on  the  caps  for  p-directed  currents. 

-  Array  containing  the  values  of  the  radially  directed  triangle 
functions  on  the  caps  for  t-directed  currents. 

-  Input  array  indicates  t-excitation  on  slot  antenna  (page  14) . 

-  9  angle  for  the  incident  wave  (degrees).  THI  «  6^  in  Equation 
(43). 

-  9  angle  for  the  radiated  fields  (degrees).  THR(K)  •»  in 
Equation  (32). 

-  0  angle  for  the  scattered  fields  (degrees)  [Equation  (43)]. 

-  Array  containing  the  values  of  the  junction  half-triangle 
functions. 

-  Array  containing  the  values  of  TC.  The  storage  method  is  the  same 

P 

as  for  T. 

-  Array  containing  derivatives  of  the  half-triangle  functions  in  the 
TBE  array  with  respect  to  Z. 

-  Array  containing  derivatives  of  the  half-triangle  functions  in  the 
TCE  array  with  respect  to  p. 

-  Array  containing  derivatives  of  the  triangle  functions  in  the  TCR 
array  with  respect  to  p. 

-  Array  containing  derivatives  of  the  triangle  functions  in  the  TCT 
array  with  respect  to  p . 

-  Array  containing  derivatives  of  the  half-trianglu  functions  in  the 
TJ  array  with  respect  to  wire  length. 

-  Array  containing  derivatives  of  the  triangle  functions  in  the  TW 
array  with  respect  to  wire  length. 

-  Array  containing  the  values  of  the  triangle  functions  on  the  wires. 

-  Array  of  weights  to  be  used  for  t-dlrected  impedance  loads  on 
either  the  BOT  or  caps. 
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TZ  -  Array  containing  the  values  of  T  .  The  storage  method  is  the  name 

V 

as  for  T. 

U  -  Imaginary  number  i. 

UXJ,U¥J,JZJ  -  Input  arrays  containing  the  x,  y,  and  z  components, 

respectively,  of  the  normal  vector  to  each  junction  disk. 

UXJ 1 , UYJ 1 ,UZJ 1  -  Arrays  containing  the  x,  y,  and  z  components,  respectively, 

of  one  of  the  orthonormal  vectors  on  each  junction  disk. 

UXJ2,UYJ2,UZJ2  -  Arrays  containing  the  x,  y,  and  z  components,  respectively, 

of  ore  of  the  orthonormal  vectors  on  each  junction  disk. 

VB  -  Array  of  BOT  voltages  corresponding  to  mode  N.  VB(K)  contains  t- 

directed  voltages  VC.  on  triangle  function  K.  VB(K  +  NM)  contains 
ni 

z-directed  voltages  on  triangle  function  K  [Equation  (37)]. 

W  -  Array  of  wire  voltages. 

XB  -  Input  array  of  x  coordinates  for  BOT  (page  9). 

X31(I)  -  x  coordinate  for  generating  curve  segment  I  (meters). 

XC  -  x  coordinate  of  cap  center  (meters). 

XJ( I)  -  x  coordinate  of  the  I-th  junction  point  (meters). 

XTEST  -  Input  variable  for  near-field  test  point  x'  in  r'  (page  16). 

XW  -  Array  of  x  coordinates  for  the  wires  (meters). 

XWl(I)  -  x  coordinate  of  wire  segment  I  (meters). 

Y  -  Array  containing  the  Ym  n  submatrix  [Equation  (41)].  In  the 

near-field  analysis,  however,  Y  contains  the  measurement  matrix  ZM. 
Ym  n  is  stored  by  columns. 

YB  -  Input  array  of  y  coordinates  for  BOT  (page  9). 

YBl(I)  -  y  coordinate  for  generating  curve  segment  I  (meters). 

YC  -  y  coordinate  of  the  cap  center  (meters). 

YJ(  .)  -  y  coordinate  of  trie  I-th  junction  point  (meters). 

YTECT  -  Input  variable  for  near-field  test  point  y'  in  r'  (page  16). 

YW  -  Array  of  y  coordinates  for  the  wires  (meters). 

YW1(I)  -  y  coordinate  of  wire  segment  1  (meters). 

Z  -  Array  containing  the  impedance  matrix  stored  by  columns. 

ZBE  -  Array  of  z  coordinates  for  the  edge  terra  on  the  BOT  surfaca 

(meters) . 

ZC(I)  -  z  coordinate  of  cap  I  (meters). 

ZE(I)  -  z  coordinate  where  the  edge  term  for  cap  I  terminates  on  the 

BOT  surface  (meters). 
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-  Input  array  Indicates  z-excitation  on  slot  antenna  (page  14). 

-  z  coordinate  of  the  I-tb  junction  point  (meters). 

-  Array  containing  the  electric  and  magnetic  modal  current  coeffi¬ 
cients  for  near-tield  calculation  [Equations  (47),  (52)].  ZK  useB 
the  same  storage  location  as  the  Y  matrix,  which  is  of  o.'der  LS, 
stored  by  columns.  Rows  1  through  3  of  ZM  contain  the  K-th  modal 
current  coefficients  for  the  electric  near-field  components  in  the 
x,  y,  and  z,  respectively,  at  the  point  (XTEST,  YTEST,  ZTEST). 
Similarly,  rows  4  through  6  of  ZM  contain  the  M-th  model  current 
coefficients  for  the  magnetic  near-field  components  in  the  x,  y, 
and  z  directions,  respectively. 

-  Input  variable  for  near-field  test  point  z'  in  r'  (page  16). 

-  Array  of  z  coordinates  for  the  wires  (meters). 

-  z  coordinate  of  wire  begment  I  (meters). 

-  Array  of  weights  to  be  used  for  z-directed  impedance  loads  on  the 


ZO'I) 

Z1(I) 


BOT. 

-  Starting  z  coordinate  for  slot  antenna  I  (meters). 

-  Ending  z  coordinate  for  slot  antenna  I  (meters). 
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